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ABSTRACT 

We examine the relationship between Type la Supernova (SN la) Hubble residuals and the properties 
of their host galaxies using a sample of 115 SNe la from the Nearby Supernova Factory (SNfactory). 
We use host galaxy stellar masses and specific star-formation rates fitted from photometry for all 
hosts, as well as gas-phase metallicities for a subset of 69 star-forming (non-AGN) hosts, to show that 
the SN la Hubble residuals correlate with each of these host properties. With these data we find 
new evidence for a correlation between SN la intrinsic color and host metallicity. When we combine 
our data with those of other published SN la surveys, we find the difference between mean SN la 
brightnesses in low and high mass hosts is 0.077 ± 0.014 mag. When viewed in narrow (0.2 dex) bins 
of host stellar mass, the data reveal apparent plateaus of Hubble residuals at high and low host masses 
with a rapid transition over a short mass range (9.8 < log(M*/M Q ) < 10.4). Although metallicity has 
been a favored interpretation for the origin of the Hubble residual trend with host mass, we illustrate 
how dust in star-forming galaxies and mean SN la progenitor age both evolve along the galaxy mass 
sequence, thereby presenting equally viable explanations for some or all of the observed SN la host 
bias. 
Subject headings: supernovae: general — cosmology: dark energy 



1. INTRODUCTION 

Type la Supernovae (SNe la) are excellent tools for 
measuring cosmological distances because of their bright - 

1 Physics Division, Lawrence Berkeley National Laboratory, 
1 Cyclotron Road, Berkeley, CA, 94720. 

2 Department of Physics, University of California Berkeley, 
366 LeConte Hall MC 7300, Berkeley, CA, 94720-7300. 

3 Research School of Astronomy and Astrophysics, The Aus- 
tralian National University, Mount Stromlo Observatory, Cotter 
Road, Weston Creek ACT 2611 Australia. 

4 ARC Centre of Excellence for All-sky Astrophysics (CAAS- 
TRO). 

5 Laboratoirc de Physique Nucleaire et des Hautes Energies, 
Universite Pierre et Marie Curie Paris 6, Universite Paris Diderot 
Paris 7, CNRS-IN2P3, 4 place Jussieu, 75252 Paris Cedex 05, 
France. 

6 Department of Physics, Yale University, New Haven, CT, 
06250-8121. 

7 Physikalisches Institut, Universitat Bonn, Nufiallee 12, 
53115 Bonn, Germany. 

8 Universite de Lyon, F-69622, Lyon, France ; Universite 
de Lyon 1, Villeurbanne ; CNRS/IN2P3, Institut de Physique 
Nucleaire de Lyon. 

9 Computational Cosmology Center, Computational Research 
Division, Lawrence Berkeley National Laboratory, 1 Cyclotron 
Road MS 50B-4206, Berkeley, CA, 94720. 

10 Centre de Recherche Astronomique de Lyon, Universite 
Lyon 1, 9 Avenue Charles Andre, 69561 Saint Genis Laval Cedex, 
France. 

11 Centre de Physique des Particules de Marseille , Aix- 
Marseille Universite , CNRS/IN2P3, 163, avenue de Luminy - 
Case 902 - 13288 Marseille Cedex 09, France. 

12 Tsinghua Center for Astrophysics, Tsinghua University, 
Beijing 100084, China. 

'■* Center for Cosmology and Particle Physics, New York Uni- 
versity, 4 Washington Place, New York, NY 10003, USA. 

14 National Astronomical Observatories, Chinese Academy of 
Sciences, Beijing 100012, China. 
Deceased. 



ness {Mb ~ —19.5) and low intrinsic luminosity disper- 
sion (~ 0.35 mag). SNe la exhibit even tighter dispersion 
(~ 0.15 — 0.18 mag) after application of empirical cor 
rection for the shape of their light curv es (Phillips 1993 
Hamuy et a l.||T995l [19961 IRiess et al"l|l9951 IPerlmutter 
et al.|1997p and their color ( (Riess et al. 1 1996} |Tripp|l 998[ ) 
These qualities make them invaluable t ools for m easuring 
the expansion history o f the Universe ( Riess et aL||1998 



r- 



et 



Perlmutter et al. 19991, providing the tightest possible 
constraints on the matter and dark energy densities as 
well as the dark energy equation of state parameter w 
'e^ 



1998; Perlmutter et al. 1999; Garnavich 
. .. aL||1998i |Knop et al.||lMWl IRiess et al.||2004| \'M 
Sullivan et al.||2011||Suzuki et al.||2012| ). 



Riess et al. 
1998 



r 



Despite the successful empirical utility of SNe la, a 
thorough census of their stellar progenitors has yet to 
be established. However, some recent studies have be- 
gun to shed light on the likely progenitors of a few in- 
dividual SNe la . Pre-explosio n Hubble Space Telescope 
(HST) imaging ( |Li et al.|20fl|) of the v ery nearby SN la 
SN 2011fe in M101 fllMugent et al.|2Ull| ) ruled out all red 
giant companions to the exploding star, while the first 
photometric measurement of SN 2011fe only 4 hours af- 
ter explosion ( |Bloom et al.|2012| confidently constrained 
the exploding star to be a white dwarf and forbade most 
main sequence companions. HST i maging of a SN la rem 



nant in the Large Magellanic Cloud ( Schaefer & Pagnotta 



2012 1 ruled out any possible non-degenerate companion 
to the exploding star (i.e. th e single degenerate, SD, see 



nario 



Whelan & Iben 1973), thereby strongly favoring 



a progenitor system comprising two white dwarfs (th e 
double degenerate, DD, scenario Iben k. Tutuk ov 1984). 
Conversely, some SNe la such as SIN 2002ic (| Wood- Vaseyj 



et al.|2004||Hamuy et al. 12003) ), SN 2005gj (Alderingetal 
2006f l^ljeto et al.||2007| )7"or recently PTFllkx ( |Dilday 
et al.||2012| ) have shown clear signs of interaction with 
some circum-stellar material, strongly favoring a SD pro- 
genitor scenario for these SNe. 

Indeed while these individual cases are very exciting 
and mark a significant step forward in the study of SN la 
progenitors, it is important to build a statistically large 
sample of SN la progenitor detections in order to study 
the influence of properties such as stellar age and metal- 
licity on the luminosities and light curve shapes of the 
resultant SNe la. Thus while interesting these individ- 
ual cases do not yet present a complete picture of SN la 
progenitors, and concerns remain that the lower metallic- 
ity and younger ages of high-redshift SN la progenitors 
could potentially alter their luminosities in a way that 
could bias cosmological parameters- 
Se veral recent studies (|Kelly et al.|2010 
iTollLampeitl et al.|2010||Gupta et al.|20 



Sullivan et al. 
'!) observed 



2010[|Lampeitlet al.|2010[|Gupta 
that the corrected brightnesses of SNe la correlated with 
the stellar masses of their host galaxies, such that SNe la 
in high-mass host galaxies are brighter than SNe la in 
low-mass galaxies after application of the usual stretch- 
and color-based b rightness correction techniques. |Gal- 



et al. 



lagher et al. (2005) found tentative evidence for a similar 
trend with respect to host galaxy gas-phase metallicity 
in star-forming hosts which was recently confir med with 
higher signi fican ce in the SDSS-SN sa mple by |Konishi 



2011|) and|D'Andrea et al.| ( 2011[ ) (though see |Gal 
lagher et al.|2008| who find, tentative evidence for an op- 
posite trend in passive hosts), and for photometrically- 
inferred metallicities by ?. The origin of this trend is 
of paramount concern for SN la cosmologists. Galax; 



iaiaxy 
mass correlates with metallicity ( Tremonti et al. [20041) 
and stellar age ( Gallazzi et al.|2005p , posing the possibil 
ity that this trend could be driven by SN la progcnito: 



ity that this trend could be driven by SN la progenitor 
metallicity or age. The average stellar ages and metal- 
licities of galaxies evolve with redshift, implying that the 
average corrected SN la brightnesses at higher redshift 
will be fainter than in the local universe if the observed 
host bias is indeed driven by progenitor age or metallicity. 
Though this trend is not strong enough to negate the evi- 
dence for cosmic acceleration, it could potentially bias es- 
timation of cosmological parameters, especially the dark 
energy equation of state parameter w and its evolution 
with redshift. 

Coincident with the concerns about biased cosmologi- 
cal results is the desire to find a means to correc t for this 
observed SN la host bias. Some authors (e.g. Sullivan 



et al. 2010 Lampeitl et al. 2010) have suggested using 
host galaxy mass as a third SN la brightness correction 
parameter (after stretch and color). This tactic might 
ameliorate any progenitor-driven luminosity effects but 
carries some important complexities to track, such as 
galaxy stellar population systematics and evolution of 
galaxy scaling relations with redshift (see Section [Fj] for 
details). We show later that these considerations are 
likely to be subdominant in the current SN la error bud- 
get, and are small compared to the magnitude of the bias 
the host term is intended to correct, but will add to the 
complexity of tracking systematic errors in SN la cos- 
mological analysis. This further motivates the search for 
the source of the observed host bias in order to poten- 
tially identify a correction technique that minimizes the 



associated host galaxy systematics. 

In this work we use SN la light curve and host galaxy 
data from the Near by Supernova Factory (SNfactory - 
Aldering et al.|2002[ ) to add to the body of observational 
data supporting the existence of a host bias. We then 
examine possible origins of the trend and prospects for 
using host-based SN la luminosity corrections. This pa- 
per is part of planned series of papers examining the 
host galaxies of SNe la from the SNfactory. The main 
companion to this paper is Childress et al. 2013a (here- 
after Paper I) which presents the main host galaxy data 
and derivation of galaxy physical properties (mass, star- 
formation rate, metallicity) from those data. Paper III 
(Childress et al. 2013c, in prep.) explores the implica- 
tions of SN la intrinsic color variability and reddening by 
host galaxy dust using Monte-Carlo simulation. Finally, 
Paper IV (Childress et al. 2013d, in prep.) fits the dis- 
tribution of SN la host galaxy masses by coupling SN la 
delay time distributions to known properties of galaxies 
in the local Universe. 

This paper is organized as follows. In Section [2] we 
summarize the SN la light curve and host galaxy data 
used in this analysis. We examine trends of stretch- and 
color-corrected SN la Hubble residuals with host galaxy 
properties in Section |3j as well as Hubble residuals de- 
ri ved using the spec troscopic standardization technique 
of Bailey et al. (2009). We discuss our results in the con- 
text of multiple similar studies and calculate the struc- 
ture of the host mass Hubble residual trend in Section |U 
We present several physical processes which may drive 
this trend in Section [5) In Section [6] we discuss the im- 
plications for SN la cosmology analyses if the observed 
host bias is not corrected, as well as the systematic errors 
introduced by the application of host-based SN la lumi- 
nosity corrections. Finally, we present our conclusions in 
Section [71 

Host galaxy physical parameters are derived assuming 
a standard ACDM cosmology with fl\ — 0.7, CIm = 0.3, 
H Q = 70 km s _1 Mpc -1 (i.e. h 70 = 1), with stel- 
lar masses and star-formation rates computed using a 



Chabrier p003 ) initial mass function (IMF), as detailed 
in Paper I. Cosmological parameters for the SNfactory 
SN la sample are left intentionally blinded to preserve the 
impartiality of forthcoming SNfactory cosmology analy- 
ses. 



2. SN AND HOST GALAXY DATA 

2.1. SN la Data 

The SNe la in this sample were discovered or followed 
by the SNfactory between 2004 to 2010. Of the 119 fully 
processed SNe la with good light curves in the SNfac- 
tory sample, 98 were discovered in the SNfactory search 
conduct ed from 2004 to 200 8 using the QUEST-II CCD 
camera ( Baltay et al.||2007 l on the Samuel Oschin 1.2m 
Schmidt telescope on Mount Palomar, California. The 
SNfactory search found many SNe of all types in the low 
redshift universe, focusing on finding SNe la in the red- 
shift range 0.03 < z < 0.08. Additional SNe la presented 
here were discovered and publicly announced by other 
nearby searches such as the Li ck Observatory Supernova 
Search (Filippcnko et al.j2001[) or the Palomar Transient 



Factory ( fKau et al.||2009l 
The SN la light curve data used in this analysis are de- 



rived from flux-calibrated spectral time series obtained 
by SNfactory wit h the SuperNova Integral Field Spec- 
trograph (SNIFS, JLantz et al.|2004| on the University of 
Hawaii 2.2m telescope on Mauna Kea. SNIFS is designed 
for observations of point sources on a diffuse or struc- 
tured background, and simultaneous observation of field 
stars in multiple filters allows for spectrophotometry ob- 
servations throu gh thin clouds. Reduct ion of SNIFS data 
is described in Aldering ct al. (2006) and upda ted in 
Scalzo et all ( |2010p , |Bongard et aTTj2011[ ), and |Buton 



etal.||2012| . 

Photometry was synthesized from SNIFS spectropho- 
tometry at all epochs using non-overlapping boxcar filter 
functions whose wavelength ranges correspond roughly to 
Johnson B, V , and R. These light curves were synthe- 
sized in the observer frame then fitted using versi on 2.2 
of the SALT2 light curve fitter (|Guy et al.j[2007 



20101 
as well 



to derive the peak magnitude in B-band (tub 

as the light curve width (xi) and color (c) parameters 

Example light curves for interestin g SNe la from the 



SNfac tory data set can be found in Scalzo et al. (2010 
2012|. Details on our light curve synthesis and fitting 
techniques can be found in Pereira et al. (2013) as ap- 
plied to SN 2011fe. 

Stretch- and color-corrected Hubble residuals and their 
errors were calculated in a blin ded cosmological a nalysis 
similar to that presented in |Bailey et al. (2009). The 



sample of SNe la ex amined here includes the sample of 



Chotard et al. (2011) but is slightly expanded to include 
some additional SINe la. For all SNe the structured host 
galaxies have been subtracted from SNIFS data using 
the method of [Bongard et al. (2011 1. Full details on the 
cosmological fitting for this sample will be presented in 
Aldering et al. (2013, in prep.). 



Later aspects of this analysis (Section 3.2) will uti- 
lize Hubble residuals derived from SN la brightnesses 
corrected using the spectroscopic stand ardization tech- 



nique developed by Bailey et al. (2009). In that work 



we showed that the spectral flux ratio 7^642/443; defined 
as the ratio of the SN la flux at 642 nm divided by the 
flux at 443 nm (with 2000 km/s binning) , correlated very 
strongly with the raw SN la Hubble residuals (i.e. with- 
out stretch and color corrections). We repeat the cal- 
culation of 7?. 6 42/443 for the subset of this SN la sample 
with observations within 2.5 days of maximum light. We 
then calculate a (blinded) luminosity correction factor 7 
and apply the luminosity correction factor 77^642/443 to 
the observed peak brightnesses of our SNe la. 

2.2. Host Galaxy Data 

Host galaxy stellar masses, star formation rates (SFRs) 
and gas-phase metallicities are derived in Paper I. We 
here briefly summarize the methods presented there and 
point to that work and references therein for further de- 
tails. 

Host galaxy stellar masses and star formation rates 
(SFRs) were derived by comparing observed multi-band 
photometry to stellar population synthesis (SPS) models. 
Photometric data for SNfactory SN la host galaxies was 
gathered from public sources as well as our own targeted 
observations. Optical photomet ry was collected f rom the 
Sloan Digital Sky Sur vey (SDSS | York et al. | 2000| ) Eighth 
Data Release (DR8, |Aihara et al.||201l[ ). MR images 
from 2MASS dskrutskie et al. 2006) were obtained at the 



NASA/IPAC Infrared Science Archive (IRSApj) , with 
NIR data fo r some SN la hosts from UKIDSS (L"a wrence| 
et al.||2007| ). UV data were obtained from the GALEX 



online data archive at MAS Tp 7 ] For very faint hosts (typ- 
ically for irig > 19.0) with STJSS photometry and those 
hosts outside the SDSS footprint, we used our instru- 
ment SNIFS in imaging mode to obtain optical images 
in standard ugriz filters. For some ho sts, g-band pho- 
tometry was obtained with Keck LRIS ( Oke et aL]|1995 l 
prior to spectroscopic observations of the hosts, and was 
later zero-pointed to either SDSS or SNIFS photometry. 
Magnitudes for our hosts in each photometric band were 
measured by perf orming common apertu re photometry 
using SExtractor ( Bertin fc Arnoutsfl 996) in dual image 
mode, and corrected tor foreground Milk y Way redden- 
ing using the dust maps of Schlc gel et al. ( 1998 ) with the 
reddening law of |Cardelli et al.| ( |l989[ ). 

Stellar masses and SFRs were calculated from photo- 
metric data using a Bayes i an me thod very similar to that 
presented in Salim et al. (2007 and references therein). 
We generated a suite ot model galaxies with physically- 
motivated star-formation histories, reddened by dust fol- 
lowing a distribution derived empirically from the SDSS 
main galaxy sample. Central values and errors for SN la 
host galaxy mass-to-light ratios and specific SFRs (sSFR 
- the SFR per unit mass) were calculated as the median 
and ±1(T values of the respective probability distribu- 
tion functions for those quantities calculated from the 
weighted (by \ 2 matching to photometry) values of the 
models. Total stellar masses were computed from the 
calculated mass-to-light ratios and the observed galaxy 
magnitude (here in g-band). Host sSFRs used in this 
analysis are the average over the past 0.5 Gyr. Our 
method was validated by examining both the ability to 
recover input galaxy parameters for our model galaxies, 
as well as with a sample of real galaxies from SDSS with 
galaxy mass and sSFR v alues me asured by the M PA- JHU 
group (Kauffmann et al. 2003; B rinchmann et aL| [2004). 
We found very good agreement in both tests, indicating 
our galaxy photometry fitting method is both accurate 
and consistent with previous studies. 

Gas-phase metallicities, specifically the logarithmic 
gas-phase oxygen abundance 12 + log(0/H), for hosts ex- 
hibiting star-formation were derived from emission line 
fluxes in galaxy spectra. Our science goals required the 
systemic host galaxy redshift as measured at the host 
galaxy core. Though the IFU capabilities of SNIFS some- 
times resulted in galaxy signal being present in SNIFS 
data cubes, the small field of view (6"x6") required addi- 
tional longslit spectroscopy observations to successfully 
include the galaxy core of most SN la hosts. Longslit 
spectra for our SN la host galaxies were obtained dur- 
ing numerous observing runs at multiple telescopes from 
2007-2011. The instruments used were the Kast Double 



Spectrograph (Miller & Stone 1993) on the Shane 3-m 
telescope at Lick Observa tory, the Low Resolution Imag- 
ing Spectrometer (LRIS - Oke et aLJ1 995) on the Keck I 
10-m telescope on Mauna Kea, the R-G Spectrograph on 
the Blanco 4-m telescope at Cerro Tololo Inter- American 
Obser vatory, the Goodman High Throughput Spectro- 
graph ( Clemens et al.|2004 ) on the Southern Astrophys- 



16 http://irsa.ipac.caltech.edu 

17 http://galex.stsci.edu 



ical Research (SOAR) 4-m telesc ope on Cerro Paction, 
and GMOS-S ( pavieset aL]ll997[ ) on the Gemini-S 8-m 
telescope on Cerro Paction. Additio nally, some hosts ha d 
spectra available from SDSS DR8 (lAihara et al. 2011). 
In this work we use only longslit and SDSS host spec- 
tra, with all longslit spectra extracted from apertures 
dominated by light in the host core. While our host 
spectroscopy sample may have slightly different aper- 
tures, our focus on the galaxy core for all spectra en- 
ables a galaxy metallicity measurement consistent with 
large surveys such as SDSS. Details on instrument con- 
figurations, data reduction, and individual observations 
are presented in Paper I. 

Emission line fluxes were measured by fitting each 
galaxy spectrum as a series of Gaussian emission line pro- 
files superposed on a linear combin ation of simple stel- 



lar po pulation (SSP) spectra from [Bruzual fc Chariot 
(2003) spaced uniformly in logarithmic stellar age, loT 



lowing the same procedure as that of |Tremonti et al.| 
(2004). Emission line fluxes were corrected for inter- 
nal reddening within the h ost galaxy by employing th e 
Balmer decrement method ( Osterbrock fc Ferland|2 006), 
and galaxies whose e mission line fluxes were contami 



nated by AGN activity (Baldwin et al. 1981 Kewley et al. 
20061 were cut. Initial mctallicities tor "high" metallicity 
hosts (defined by log(NH/ Ha) > —1.3) were calc ulated 
using the "N2" method of |Pettini fc Pagel] p004| , while 
for "low" metallicity hosts initial metallicitie s were cal- 
cula ted us ing the "R23" method of |Kobulnicky fc Kew- 



ley (120041) . All metal licities were then converted to the 



1'remonti et al. ( 20041) scale using the con version formu- 
lae presented in Kewley & Ellison (2008). This enables 
us to use the strongest emission lines across the full range 
of metallicities while retaining a common scale (see dis- 
cussion in Paper I). The systematic uncertainty in such 
metallicity conversions is typically about 0.06 dex, so we 
have added this value in quadrature to the metallicity 
measurement errors for our SN la hosts. 

2.3. SN la Light Curve vs. Host Properties 

We plot the SN la light curve parameters against the 
properties of their host galaxies in Figure [11 Trends of 
SN la properties with the properties of their host galaxies 
are important for detecting underlying physical mecha- 
nisms driving the shape of SN la light curves or their 
colors. 

Early studies of SN la host galaxies found qualita- 
tive evidence for a correlation between the observed peak 
magnitude (before correction for light curve width and 
color), light curve decline rate, and expansion velocity of 
an SN la with th e morphologica l type of its ho st galaxy 
(iFlhpp cnko 1989; Branch & van den Bergh|19"93J |Hamuy| 



et a 



1996 Gallag her et al.||2005p such that brighter 
slower declining SNe la preferentially occur in later type 
(spiral and irregular) galaxies, while fainter faster declin- 
ing SNe la preferentially occur in earlier type (elliptical 
and SO) galaxies. Recent analyses of large samples of 
SN la host galaxy masses have yielded the similar result 
that the observed brightnesses of SNe la c orrelate with 



the stellar masses of the ir host galaxy (e.g. Howell et al 



fai 



2009 Neill et al. 20091, such that more massive hosts 



produce preferentially fainter SNe la. 

Our sample of SNfactory SNe la shows a similar trend 
of light curve widths with the stellar masses of their host 



galaxies. We also detect trends of light curve width with 
galaxy sSFR and metallicity. Such correlations are ex- 
pected given that these quantities correlate with stellar 
mass in normal galaxies. What is perhaps most interest- 
ing about these figures is that rather than exhibiting a 
tight correlation of stretch with host mass or metallicity, 
the data instead indicate the presence of both high and 
low stretch SNe la in high mass (metallicity) hosts, while 
low mass (metallicity) galaxies seem to host exclusively 
intermediate to high stretch (i.e. brighter) SNe la. 

An examination of SN la colors in the SNfactory sam- 
ple shows that very red SNe la (e.g. c > 0.2) only occur in 
high mass (and high metallicity) hosts. This is perhaps 
to be expected as these SNe la are predicted to suffer 
from extinction by foreground dust, which is more abun- 
dant in high mass (and high metallicity) st ar-forming 
galaxies (Lee et al. 2009 Garn & Best 2010). Interest- 
ingly, we notice an apparent trend of SiN la color with 
host metallicity, even excluding the apparently highly- 
reddened outliers. We will revisit this trend and its im- 
plications in Section [5] 

3. SN la BRIGHTNESSES AND HOST GALAXY 
PROPERTIES 

Using the SN la and host galaxy data described above 
(see Section [2]), we investigate the correlation of SN la 
host galaxy properties with the brightnesses of SNe la 
after the application of the standard stretch- and color- 
based luminosity corrections. In the first part of this 
Section we will use the SALT2 Hubble residuals, then 
examine the Hubble residuals derived usi ng the spec- 
trosco pic luminosity correction method of |Bailey et al. 
(2009) and their correlation with host properties. 

Before beginning the analysis, we briefly discuss our 
procedure for weighting the data. We require that fitted 
models have reasonable goodness-of-fit, as judged using 
the value of x 2 P er degree of freedom, xt- F° r most su- 
pernova Hubble residual analyses it has been necessary 
to add an extra intr insic dispersion to obtain \t ~ 1 ( e -g- 
Conley et aLl|2011[ ). The SALT2 version 2.2 lightcurve 
model has an associated uncertainty that is propagated 
into the uncertainty on the Hubble residuals. For the 
SNfactory dataset the intrinsic uncertainty on the SALT2 
color dominates the error budget, and leads to very sim- 
ilar uncertainties for all of our SNe. Ultimately, these 
uncertainties give xt ~ 0-94 (for 116 degrees of free- 
dom), which is acceptably close to 1. The size of the 
observed SN la host mass bias is small compared to the 
observed scatter of SN la Hubble residuals and thus has 
little affect on the goodness-of-fit. In light of these con- 
siderations, for the analysis in this Section we do not add 
an intrinsic dispersion to our Hubble residuals when fit- 
ting linear models. However, for binned means we derive 
uncertainties based on the reduced weighted RMS (the 
weighted RMS divided by the square root of the number 
of degrees of freedom) for the data in each bin, which 
effectively accounts for the observed scatter of the data 
even if the uncertainties assigned to each SN were to be 
slightly misestimated. 

3.1. Hubble Residuals vs. Galaxy Properties 

In this analysis we use the Hubble residuals for 119 
SNe la from SNfactory after corrections have been made 
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Figure 1. Light curve width (SALT2 x\ - top row) and color (SALT2 c - bottom row) for our SNe la as compared to the properties of 
their hosts: stellar mass (left column), specific star formation rate (middle column), and gas phase metallicity (right column). 



for light curve width (SALT2 xi) and color (SALT2 c). 
Of these 119 SNe la, 115 have host stellar mass and sSFR 
estimates from photometry, 2 are apparently hostless, 
and 2 are lacking sufficient host photometry data. Of the 
119 SNe la considered here, a total of 69 (all of which 
are in the subset with host masses) have good gas-phase 
metallicities (i.e. are star forming and have emission line 
fluxes consistent with star formation rather than AGN 
activity - see Section [2] for details) . 

In Figure[2]we plot the SALT2 Hubble residuals against 
host galaxy stellar mass, sSFR, and metallicity respec- 
tively. For visual aid we show the bin-averaged values of 
the SN la Hubble residuals in bins of host mass (sSFR, 
metallicity), as well as the best fit linear trend and mean 
residuals split by host mass (sSFR, metallicity). 

Our data indicate a correlation of SN la Hubble resid- 
uals with all three host galaxy properties: mass, sSFR, 
and metallicity. We quantify this observed trend with 
two metrics: a linear trend in Hubble residuals vs. host 
stellar mass (sSFR, metallicity), and the difference be- 
tween the average Hubble residuals when SNe la are 
split into two bins corresponding to high- and low-mass 
(sSFR, metallicity) hosts. The linear trend is fit by min- 
imizing the sum of \ 2 m the two-dimensional Hubble- 
residual host mass (sSFR, metallicity) plane via an or- 
thogonal distance regression. We note the slopes and 
X 2 values for this method are extremely close to a sim- 
ple one-dimensional linear fit obtained by minimizing the 
Hubble residual \ 2 1 due to the very shallow slope of the 
trend. For this and following sections the uncertainty 
quoted on the "step" value corresponds to the quadra- 
ture sum of the uncertainties on the weighted Hubble 
residual mean values for each side of the step. We con- 
firmed our confidence intervals on the step by randomly 



re-assigning host masses to Hubble residual values, and 
examining the distribution of step values from many re- 
alizations of this process. 

We summarize the results of these fits in Table [T] In 
this table we also present the \ 2 values for the linear fits 
and split samples (which have two degrees of freedom less 
than the number of SNe for each subset) . For reference, 
the x 2 °f the Hubble residuals without any linear trend 
or step is 109.4 for the full sample of 115 SNe la, and 69.5 
for the 69 SNe la with host metallicities. The linear fit 
of Hubble residuals against host mass has a significant 
slope. Likewise, the difference in Hubble residuals be- 
tween high-and low-mass-hosted SNe la (whi ch we chose 
to sp lit at log(M*/M Q ) = 10.0 following Sullivan et al. 
2010 1 shows a significant step in corrected SIN la bright- 
nesses. We find that low- and high-mass-hosted SNe la 
have brightnesses that differ by 0.085 ± 0.028 mag after 
stretch and color corrections have been applied, such that 
SNe la in high mass hosts are brighter after correction 
than those in low mass hosts. Our data also indicate that 
corrected SN la brightnesses differ for SNe la in hosts 
w ith different metallicit y (split at 12 + log(0/H) = 8.8 as 
in D'Andrea et al.|2011 ) . Linear or step- wise trends with 
st ar- format ion int ensity (split at log(sSFR) = -10.3 as 



Sullivan et al. 2010[ ) are weaker, although consistent 



with previous studies. These results imply that applying 
the canonical SN la standardization techniques to SNe la 
at all redshifts will result in biased cosmological parame- 
ters as the SN la environments evolve in metallicity, and 
perhaps star- format ion activity, as a function of redshift. 
We quantify this potential bias in Section [6j 

To demonstrate the robustness of the observed host 
bias, we plot in Figure [3] the SN la Hubble residuals for 
our sample against two model-free observational quan- 



Table 1 

Hubble Residual Trends with Host Properties 



Host 
Property 



Residual N$Ne 

Type 



Linear Trend Linear 



K 



(mag/dex) 



Split N lo /N hi 
Value 



Hubble Residual Step HR Step 



Step (mag) 



X 



Significance 



Mass 

sSFR 

Metallicity 



Stretch+Color 115 
Stretch+Color 115 
Stretch+Color 69 



-0.043 ±0.014 

0.031 ± 0.008 

-0.106 ±0.043 



100.5 

104.6 

62.1 



10.0 52 / 63 

-10.3 46 / 69 

8.8 30 / 39 



0.085 

-0.050 

0.103 



± 0.028 98.8 
±0.029 105.0 
± 0.036 60.0 



3.0c- 
1.7ct 
2.9ct 



M a 


Stretch+Color 


115 


0.021 ± 0.008 


101.5 


-19.8 


57 / 58 


-0.071 ± 0.028 


101.6 


2.5ct 


9-i 


Stretch+Color 


112 


-0.059 + 0.030 


103.0 


0.82 


56 / 56 


0.069 + 0.029 


100.9 


2.4ct 


Mass 


Flux Ratio 


94 


-0.016 + 0.014 


92.4 


10.0 


38 / 56 


0.050 + 0.028 


90.6 


1.8a- 


Mass 


Stretch+Color 


94 


-0.042 ± 0.016 


77.7 


10.0 


38 / 56 


0.066 ± 0.030 


78.8 


2.2ct 


sSFR 


Flux Ratio 


94 


0.047 + 0.015 


85.4 


10.0 


39 / 55 


-0.052 ± 0.028 


87.6 


1.9c 
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Figure 2. Top: SALT2 Hubble residuals for SNfactory SNe la 
plotted versus host galaxy stellar mass (grey points). The blue 
line represents the best fit linear trend, the green points represent 
binned averages, and the thick red lines represent the averages for 
Hubble residuals split into high and low mass bins. Middle: Same 
as top, but for host sSFR. Bottom: Same as top, but for host 
gas-phase metallicity. 



tities: the host galaxy absolute magnitude (in g-band), 
and the host galaxy color (g — i). Trends in SN la Hub- 
ble residuals with these properties are again evident, and 
we report the linear trend and magnitude steps (split at 
the median value of these quantities) for these host prop- 
erties in Table [I] Note that three SN la hosts did not 
have an i-band magnitude and thus do not appear in the 
sample with host g — i colors. 

We note that these quantities are in the observer frame 
of reference and have not been K-corrected for the host 
galaxy redshift. Since this is a supernova-selected sample 
of galaxies, redshifts are uncorrelated with galaxy abso- 
lute magnitude or color and thus will not significantly 
alter the morphology of these plots. Furthermore, be- 
cause of the low redshift range of our sample (z < 0.1), 
K-corrections were always less than 0.05 mag, compara- 
ble to the size of the errorbars in Figure [31 Addition- 
ally, we have performed K-corrections for galaxy magni- 
tudes (which are model- dependent) and confirmed that 
the step in SN la Hubble residuals persists when using 
K-corrected photometry to derive host absolute magni- 
tudes and colors. 

Thus we claim that the observed bias of SN la Hubble 
residuals with the properties of their host galaxies is not 
an artifact of erroneous modeling of host galaxy physical 
properties (indeed one would not expect such errors to 
correlate with supernova Hubble residuals in such a way 
as to produce a trend), but instead reflects a real physi- 
cal relationship between the properties of SN la Hubble 
residuals with the nature of their environments. 

We now discuss briefly why the observed trend of the 
corrected SN la Hubble residuals with host mass is very 
unlikely to be due to selection bias. One possible con- 
cern would be that in high mass galaxies SN la searches 
are only finding comparatively brighter SNe la due to 
the high background brightness from the massive host 
galaxy. This is particularly concerning for high redshift 
surveys where the angular diameter of SN la host galax- 
ies is comparable to the ground-based seeing element. 
When we consider the distributions of both host galaxy 
masses and SN la stretches (i.e. x\) as a function of 
redshift for our SNfactory sample, we find no detectable 
change in the shape of these distributions as a function 
of redshift up to at least z = 0.12. This is explored in 
more detail in Childress et al. (2012c, in prep.) where we 
examine the SN la host galaxy mass distribution. Here 
we note that this result implies that SNe la from SNfac- 
tory are not preferentially missed in high mass galaxies 
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Figure 3. Same as Figure [2| but with Hubble residuals plotted 
against model-free observational data: host absolute magnitude (in 
g-band, top panel), and host color (g — i, bottom panel). 

as the redshift increases and the host background thereby 
increases due to decreased angular size. 

3.2. Spectroscopically- Corrected SN la Brightnesses 
and Host Properties 

We now revisit the above analyses with an alternate 
SN la lumino s ity sta ndardization technique developed by 
Bailey et al. ((2009), which employs spectral flux ratios 
to standardize SNe la. The objective of this Section is 
to examine whether the host bias in corrected Hubble 
residuals persists when using a spectroscopic brightness 
correction technique, or whether it is only photometric 
standardization methods that suffer this bias. 

The SN la data set for this analysis consists of the 
subset of the 119 SNe la from the previous sam ple for 
which we c an apply the correction method of iBailcy 
et al. (2009). This method requires a spectrophotomet- 
ric observation of the SN within ±2.5 days of -B-band 
maximum light (as estimated from the full light curve). 
This requirement brings the parent sample of SNe la 
down to 98, of which we have host stellar masses for 
94. The measurement uncertainties using this technique 
are small, and the resulting scatter about the Hubble dia- 
gram revealed the presence of an intrinsic dispersion com- 
ponent. Thus a blinded intrinsic dispersion was added in 



Figure 4. Same as Figure [2] but with Hubble residuals ob- 
tained using spectral flux ratio Drightncss corrections following the 
method described by [Bailey et al,|(|2009|l. 



quadrature to the measurement errors in order to pro- 
duce xt = 1 for the Hubble diagram. 

For simplicity in this section, we inspect only the 
trends of Hubble residuals with stellar mass and sSFR, 
since the sample attrition for a metallicity analysis is 
rather significant. In Figure [4 we show the flux-ratio- 
corrected SN la Hubble residua s for our 94 SNe la plot- 
ted against the stellar masses of their host galaxies. Sim- 
ilarly to Figure [2| we plot the binned average Hubble 
residuals, best fitlinear trend, and average Hubble resid- 
ual when splitting the sample by host mass. 

The results of our analysis show a trend with host mass 
similar to that found for stretch- and color-corrected 
Hubble residuals. The high- to low-mass magnitude step 
after flux-ratio correction is 0.050 ± 0.028 mag, versus 
0.085±0.028 mag from stretch- and color-corrected SN la 
brightnesses. The linear best-fit trend of flux ratio Hub- 
ble residuals with host mass is —0.016 ± 0.014 mag/dex, 
compared to —0.043 ± 0.014 mag/dex for stretch- and 
color-based Hubble residuals for the same sample. 

While at first glance this appears to be a decrease in 
the observed host bias, the improvement is less marked 
when comparing the bias on common SN la subsamples. 
Since those SNe la without a spectroscopic observation 
near maximum light were excluded from this sample, we 
repeat the original calculation of the previous section 
for stretch- and color-corrected Hubble residuals for the 
same sample of SNe la used in the flux-ratio analysis. 
Doing so we find the luminosity step for this sample is 
0.066 ± 0.030 mag, and the slope of the linear trend fit is 
—0.042 ± 0.016 mag/dex. These results are summarized 
in Table [T] along with the rest of our Hubble residual 
trend results. 

The result is similar when considering trends in sSFR. 
The flux-ratio Hubble residual step in sSFR for this sam- 
ple is -0.052 ± 0.028 mag versus -0.038 ± 0.030 mag 
for the same sample with stretch- and color-corrected 
Hubble residuals. Similarly, the slope of flux-ratio Hub- 
ble residuals with sSFR is 0.047 ± 0.015 mag/dex versus 
0.021±0.017 mag/dex for the same sample of SNe la with 
SALT2 Hubble residuals. Thus the flux-ratio-corrected 
Hubble residuals show a host bias which is consistent 



with that found using standard stretch and color correc- 
tions to SN la luminosities. 

4. SN la HUBBLE RESIDUAL HOST BIAS 
RESULTS FROM MULTIPLE DATA SETS 

Our work reinforces the results of recent SN la host 
studies which all have found that the stretch- and color- 
corrected luminosities of SNe la are biased with respect 
to the properties of their host galaxies. Such a trend 
is quite distressing for future SN la cosmology missions 
which will hunt for SNe la at very high redshifts where 
the mean stellar age and metallicity are quite different 
from their values in the local universe. The discovery of 
this SN la host bias has motivated two critical questions: 
what is the true cause of this bias, and what should be 
done to correct it? We will address each of these ques- 
tions in turn in Sections [5] and [6j but first we analyze the 
impact of our SNfactory data on their answers. In this 
Section we begin by comparing the results of our anal- 
ysis with those of previous studies, then we inspect the 
strength of the observed SN la host bias when combining 
data from multiple surveys. 

4.1. Comparison To Previous Studies 

In this section we summarize the results presented here 
and compare those to the findings of previous studies. 
We compare our resul ts to previous stud i es with respect 
to host stellar mass ([Kelly et ah] [20101 |Sullivan et al. 



2010 Lampeitl et al.|20lU[|CJupta et al.|201ip, then com- 



pare our tre nd w ith host metallicity to" t hose oflKonishi 
et al.| ( |2011[ ) and |D'Andrea et al.| ( |2011| ). In Table |] we 

summarize the results of these previous studies in terms 
of the Hubble residual slope and step values, as well as 
their consistency with the results derived in this work for 
SNfactory. For reference we also show the results for the 
combined sample analysis outlined below in Section [4. 2 1 
Two items in this table are worth noting here. The 
first is that the Hu bble residual step (with respect to 
host mass) quoted in |Kelly et al.| ( [2010[ ) comes from split- 
ting that sample at iog{M if / Mq) = 10.8, rather than at 
10.0 as was done for all other sample, and so is not di- 
rectly comparable to our SNfactory result. The second 
item of note regards the uncertainty in Hubble residual 



step (w ith respect to metallicity) for the D 'Andrea et al. 
(2011) results. These authors reported a step value of 
I\m C orr = 0.114 ± 0.023 mag, but the error bar on this 
result is derived from Hubble residuals whose uncertain- 
ties were not adjusted to give an acceptable goodness of 
fit. This choice was an attempt to avoid deweighting the 
data in the event most of the remaining scatter arose 
from host-dependent bias. However, the linear and step 
models re move little of the excess scatter in the data pre- 
sented by D'Andrea et al. (2011) so their quoted uncer- 
tainty underestimates the true scatter in the data. Given 
their sample size, we estimate th at an intrinsic scatter 



term of (Tint = 0.14 mag (cited by D'Andrea et al 



2011 



as the scatter in their Hubble residuals) would make the 
uncertainty in their result increase to 0.053 mag. 

We see from Table [2] that the SNfactory results are 
fully consistent with previous studies of the SN la host 
bias. The only ma jor exception is the disagreement of the 
Kelly et al. (20101 Hubble residuals sl ope with both our 
results and those of other authors. The |Kelly et al. ( 2010 1 



sample covered a shorter mass range than the other three 
samples, and we will show in the following section that 
their result is consistent with detailed structure that we 
have uncovered in the trend of SN la Hubble residuals 
over their shorter host mass range. 

4.2. Hubble Residual Host Trends from Combined 
Data Sets 

Our analysis of SNfactory SN la Hubble residual trends 
with host galaxy properties adds to the continually grow- 
ing body of data that indicates current SN la standard- 
ization methods produce a bias of Hubble residuals with 
respect to host galaxy properties. We now examine the 
observed bias when multiple data sets are combined, in 
part to address any concern that such a bias could be 
the result of spurious effects or inaccurate galaxy stellar 
population synthesis methods. Here we will show that 
the combination of multiple data sets reveals detailed 
structure of the Hubble residual trend with host mass. 

For this purpose, SN la Hubble residuals and host 
galaxy data have been collected from two additional un- 
targeted S N la searches: the Supernova Legacy Su rvey 



(SNLS - |Sullivan et aT] [20101 |Guy et al. [2010] ) and 



the SPSS-SIN survey QGupta~eT atpm] ), as weTTas the 



compilati on of low-redshift literature SNe la and their 
hosts of Kelly et al. (2010). We a pply a common flat 
ACDM cosmology w ithHIF = 0.271 (Sullivan et al.|2011' 



Suzuki et al. 2012 ) and SN la stand ardization parame 
ters a = 0.139 and f3 = 3.1 (Guy et al.||2010| ). We note 
that these SN la luminosity correction parameters are 
appropriate for the data here because they are derived 
for the same version (2.2) of the SALT light curve fitter 
used to fit all the data, and were derived without fitting 
simultaneously for a host mass correction term. 

Host galaxy ma s ses fo r the SNLS sample were derived 
by Sullivan et al. (2010), and SN la light curve param - 
eters and redshitts were presented in Guy et al. (2010). 
Hubble residuals for this sample were derived directly us- 
ing the aforementioned cosmology and SN la standard- 
ization parameters. Host galaxy data and SN la Hub- 
ble residuals f or the SDSS-SN sample were presented in 
Gupta et all poTT}. We applied small corrections to 



their reported Hubble residuals and errors in order to 
conv ert from their cosmolo gy and light curve parame- 
ters ( Kessler et al.[ ;2009; Holtzman et al. 2008) to ours. 
Host masses for the literature low-redshift sample were 

)|, w hile light curves were 
(2009). We use new Union3 

light cur ve fits for 61 of the 



derived in |Kelly et al 



presented in Hicken et al 



(2010) 



(Rubin et al. 2013, in prep 



62 SNe la from the Kelly et al. (2010) sample using ver 



sion 2.2 of SALT, omitting the slowlleclmer SN 2006bz. 
These updated light curve fits were performed in order 
achieve consistency with the three other samples, as the 
updated SALT color model beha ves different ly than that 



used in the light curve fits from Kelly et al. (12010). For 
SNfactory data, we use the raw blinded Hubble residuals 
(before stretch and color correction) then apply correc- 
tions for the chose n SN la st andardization parameters. 

We note that the Sullivan et al.| ( ]2010[) and Gupta et al. 
( |2011[ ) host masses are both computed using a Kroupa 
(TOT 

ChaFrier| ( |2003[ ) IMF. Fortunately these IMFs are ex- 
tremely similar, and while they may exhibit some rela- 
tive scatter (perhaps 0.05 dex) they do not require an 
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Authors 


SN 


Host 


N S Ne 


Median 


Median 


Linear Trend 


SNfactory 


Hubble Residual 


SNfactory 




Survey 


Property 




Redshift 


M or Z 


(mag/dex) a 


Agreement 


Step (mag) a 


Agreement 


Sullivan '10 


SNLS 


M 


195 


0.64 


9.8 


-0.042 ±0.013 


0.05(7 


0.08 ± 0.02 


0.15cr 


Lampeitl '10 


SDSS 


M 


162 


0.22 


10.3 


-0.072 ±0.018 


1.27o- 


0.100 ±0.025 


0.40(7 


Gupta '11 


SDSS 


M 


206 


0.22 


10.3 


-0.057 ±0.019 


0.59ct 


0.096 ± 0.028 


0.28(7 


Kelly '10 


Low-z 


M 


62 


0.03 


10.7 


-0.15 ±0.06 


1.74ct 


0.094 ± 0.045 6 


— 


This Work 


SNf 


M 


115 


0.03 


10.1 


-0.043 ± 0.014 


— 


0.085 ± 0.028 


— 


This Work 


All 


M 


601 


0.25 


10.2 


-0.037 ±0.011 


0.34(7 


0.077 ±0.014 


0.26cr 


Konishi '11 


SDSS 


Z 


72 


0.22 


8.9 


— 


— 


0.13 ±0.06 


0.39cr 


D 'Andrea '11 


SDSS 


z 


34 


0.22 


8.8 


— 


— 


0.114 ±0.053 c 


0.16cr 


This Work 


SNf 


z 


69 


0.06 


8.9 


-0.106 ±0.043 


— 


0.103 ±0.036 


— 



Kroupa (20011 IMFs and found a mean off- 
dex (which we note is much smaller than 



"Values as originally published by stated authors. 

6 Split at log(M*/M ) = 10.8 

c Including added intrinsic scatter, see text. 

offset in mass-to-light ratio zeropoints. The host galaxy 
sample of Kelly ct al. (2010) uses masses computed with 
th ejRana'fc Basu (1992) IM F, which is the defa ult IMF 
in ZPEG (Lc Borgnc & Rocca-Volmerange 2002). From 
a sample ot approximately 10,000 galaxi es from SDSS, 
we com put ed ZPEG masses using both the |Rana &: Basu| 
(1992) and 
seTof 0.02. 

typical mass uncertainties o f 0.1 dex), which we subtract 
from the Kelly et al. (2010) masses for consistency with 
our IMF and that ot SiNLb and SDSS. Our host masses 
and those reported by the aforementioned authors were 
all computed using the same flat ACDM cosmology with 
&m = 0.3 and h — 0.7, indicating our mass scales are all 
consistent. 

The combined sample of Hubble residuals as a function 
of host galaxy mass is presented in Figure [5] The total 
number of SNe la included in this analysis now totals 601 
(219 SNLS, 206 SDSS-SN, 61 literature low-z, and 115 
SNfactory). We now wish to calculate the step value and 
best fit linear trend from this combined data set, but 
once again must consider the proper weighting of the 
data. The measurement uncertainties alone yield xt ~ 2 
for the combined sample. Thus, to obtain xl ~ 1 we add 
an intrinsic dispersion term of 0.10 mag in quadrature 
with the measurement uncertainties for all samples in 
the combined dataset, which yields a favorable x 2 of 628 
for the sample of 601 SNe la. 

With this combined data set we fit for a linear trend 
using the same two-dimensional % 2 minimization method 
employed for the SNfactory linear fits in Section [3j This 
is shown as the thick dashed orange curve in Figure [51 
and has a slope of —0.037 ± 0.011 mag/dex. We then 
calculate the step in average Hubble residuals between 
SNe la in high- and low-mass hosts when splitting the 
sample in host mass at log(M*/M0) = 10.0 (this results 
in Ni ow — 255 and N^g^ — 346 SNe la in the two bins). 
We find a combined offset in corrected SN la brightnesses 
of 0.077±0.014 mag, a 5.6c detection of host bias in SN la 
Hubble residuals. 

Concurrent with the measurement of the Hubble resid- 
ual step for the combined SN la sample, we also cal- 
culate the step for SNLS and SDSS separately to com- 
pare against previously published results. We find the 
Hubble residual step is 0.065 ± 0.020 mag for SNLS 
and 0.114 ± 0.025 mag for SDSS-SN. We note that the 



step value f or SNLS we deri v e here is smaller than that 
reported in |Sullivan et al. (2010), and we confirmed 
that this is due to the use of versi o n 2.1 of the SALT 
light curve fitter in Sullivan et al. (20101. Finally we 
checke d the Hubble residual step for the |Kelly et al.| 
(|2010| sample split at their value of log(M*/Af Q ) = 10.8 



using our chosen cosmology and found this step to be 
0.108±0.045 mag, consistent with their reported results. 

4.3. Structure in the Host Mass Hubble Residual Trend 

We now use the combined SN la host mass and Hubble 
residual sample to inspect the structure of the trend of 
SN la Hubble residuals with host properties. We bin the 
sample in narrow bins of 0.2 dex in host galaxy stellar 
mass, and plot the result in Figure [6] The result of this 
analysis indicates that rather than a clear linear trend 
or step in host mass, the binned SN la Hubble residuals 
appear consistent with a plateau at both low and high 
stellar mass separated by a narrow transition region from 
log(M*/Af ) = 9.8 to log(M*/M ) = 10.4. 

To further illustrate this result, we recalculate the 
weighted mean Hubble residual for all SNe la with hosts 
more massive than log(M,/M ) = 10.4, and find that 
the difference between this value and the weighted mean 
Hubble residual for SNe la in hosts less massive than 
log(M*/M ) = 9.8 is 0.086 ± 0.016 mag. This value is 
larger than the canonical value we calculated from our 
split of all SNe la at host mass of log(M»/M ) = 10.0. 
This indicates that the Hubble residual host bias ob- 
served in the data is not best represented by a sharp step 
function and that previous measurements of the offset be- 
tween SN la luminosities in high and low-mass hosts may 
be underestimating the true value of the offset. 

We note that this result is not as clearly evident when 
combining the corrected Hubble residuals from the four 
samples without first applying a common cosmology and 
SN la correction parameters (a and /?), nor is it as pro- 
nounced if any one of the SN la samples from untargeted 
searches (SNLS, SDSS-SN, SNf) is excluded. This result 
then is found only with the most up-to-date and compre- 
hensive sample of SN la host galaxy masses and Hubble 
residuals. 

This result also elucidates the likely cause of the strong 
trend of Hubble r esiduals with host mass reported by 
Kelly et al. (2010), which we previously noted as being 
much stronger than that reported by other studies in- 
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Figure 5. Top: Com bined Hubble residuals and host galaxy m asses for the SNfactory sampl e (this work), the SNLS survey (Sulliva n et al.| 
|2010||Guy et al.|2010} , the SDSS-SN survey jGupta et al.|2011) , and literature low-z sample | |Kelly et al |2010||Hicken et al.lZuUS) ; markeoT 
as blue circles, green diamonds, magenta triangles, and red squares, respectively. Symbol size is inversely proportional to iiubble residual 
error, and open symbols represent SNe la whose Hubble residual measurement errors exceed 0.2 magnitudes. The weighted mean Hubble 
residuals for high- and low-mass-hosted SNe la are shown as thick black bars, and the best fit linear trend is the thick dashed orange line. 
Bottom: Histograms of host mass distributions for the four samples. 



eluding our own. If most of the Hubble residual bias 
between low- and high-mass hosted SNe la is driven by a 
rapid shift in mean Hubble residuals over a short range, 
then the local Hubble residual gradient can be quite 
strong. Indeed if the entire offset of 0.086 ± 0.016 mag 
transitions over a mass range of 0.6 dex, then one would 
expect a gradient of 0.14 ± 0.03 mag/dex, which is re- 
markably close to the trend reported by |Kelly et al.| 
(2010). Most of their SN la hosts were in the mass range 
ot this transition, and thus serendipitously occurred in 
the region where the host bias is strongest. 

For completeness here, we note that the weighted mean 
and RMS values have only a mild dependence on how one 
chooses to weight individual SNe la. For this Section and 
the analysis of SNfactory data (Section [3]) weights were 



assigned according to the measurement errors alone, but 
our use of the reduced weighted RMS in assessing the 
uncertainty absorbed any unaccounted-for scatter. An 
alternative approach is to add an intrinsic dispersion in 
quadrature with the measurement uncertainties. This 
preserves the relative order of weighting of SNe la but 
decreases the difference in weights between SNe la with 
high and low measurement error. We have confirmed 
that adding an intrinsic dispersion term of 0.10 mag to 
the measurement uncertainties for all SNe la results in 
negligible (< 0.5c) change to the weighted mean and 
RMS values reported here. 

5. POSSIBLE ORIGINS OF SN la HUBBLE 
RESIDUAL HOST BIAS 
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Figure 6. Binned average Hubble residuals from Figure[5]in host 
mass bins of width 0.2 dex. The highest and lowest mass Dins are 
wider in order to have at least 20 SNe in each bin. Error bars 
for each bin are the weighted RMS divided by the square root of 
N — 1 where N is the number of points in that bin. An apparent 
region of Hubble residual transition from \og,(M* / Mq) = 9.8 to 
log(M*/M0) = 10.4 is denoted by the vertical dashed lines, while 
the mean Hubble residuals of SNe la in hosts of masses above and 
below this region are denoted by the thick red bars. 

The first major question raised by the SN la luminosity 
host bias is, what is its underlying physical cause? Is it a 
true correlation of SN la brightnesses with the properties 
of their progenitors? If so, is the trend being driven 
by progenitor age or metallicity (or possibly both), or 
perhaps some other physical property of the progenitor? 
On the other hand, could this be a deficiency in the SN la 
standardization techniques? In this Section we address 
these questions by comparing the observed SN la host 
bias to a series of physical models. Specifically we seek 
a physical model which can reproduce the fast transition 
of Hubble residuals over a short mass range discovered 
when we combined multiple data sets (see Section HI). 

One candidate driver of the host bias in SN la lumi- 
nosities is erroneous color correction resulting from the 
interplay of SN la intrinsic colo r va riation and redden- 
ing by foreground dust (Section |5.1| and Figure u\. Our 
SNfactory data indicated a possible tre nd o f SN la intrin- 
sic color with host metallicity (Section 5.2 and Figure[8), 
which we show could drive a trend of Hubble residuals 
with host properties. We then discuss the likely shift of 
average SN la progenit or ag es across the range of host 
galaxy masses (Section |5.3| and Figure [9]) and its po- 
tential for driving some of the observed nost bias. Fi- 
nally we consider the effect of neutroniza tion in SNe la 



(Timmes et al. 2003 Bravo et al. 2010) and its effect 



on metallic ity- dependent stretch corrections (Section 5.4 



and Figure 10). 



Our findings ar e th en summarized in Section |5.5[ where 
we plot in Figure [TT] the observed trend of SN la Hubble 
residuals with host galaxy mass (see Section El) against 
the trends predicted from all the models presented in this 
Section. 

5.1. Competing Effects of Intrinsic SN la Color and 

Dust 

A potential culprit that could drive the observed trend 
of SN la Hubble residuals with host properties is the cor- 



rection of SN la luminosities for color. If there is signifi- 
cant variation in the intrinsic colors of SNe la, especially 
if intrinsic color affects the SN la brightness differently 
than dust (i.e. different /3), then correction of SN la 
brightnesses via a single color law will produce erroneous 
values. The magnitude and sign of the color correction 
error for each SN la is then related to the relative pro- 
portion of dust and intrinsic color variation. This can 
produce trends with host galaxy properties due to the 
corre l ation of galaxy dust content with mass (Lee et al. 



2009 Garn & Best 2010) such that more massive star 



forming galaxies are more dusty. The magnitude and 
direction of the SN la color correction error is dependent 
on the slopes of the SN and dust color laws. 

Here we illustrate this effect with a specific example 
utilizing Monte-Carlo simulations of mock SNe la. A 
key input to to simulation is the SN la intrinsic color 
distribution, for which we choose a Gaussian of w idth 
0.07 mag (e.g. |Jha et al.|[2007| |Folatelli et aTJ|2olo|) We 
further assume that the intrinsic color affects the SN la 
luminosity following some slope /3sn = 5.0 such that a 
SN la which is intrinsically Ac = 0.02 mag redder would 
be Am = /Sg^Ac = 0.10 mag fainter. 

We generate a mock sample of SNe la whose intrinsic 
colors follow this distribution and are found in galaxies 
whose masses are distributed similarly to the combined 
sample of SNe la from Section E] Utilizing the observed 
correlation of SN la stretch with host mass (see e.g. Fig- 
ure [l]), we assign a stretch to each mock SN. 

Finally we assign a host galaxy reddening to each SN la 
using the correlation of dust content with host mass in 
star-forming galaxies. The highest mass mock SN hosts 
are considered to be passive dust-free ellipticals with a 
likelihood of being passive trained on SDSS data. Then a 
foreground dust reddening of the SN was drawn from an 
exponential distribution scaled to the total host galaxy 
dust content. This critial piece accounts for the fact that 
an SN is likely to be geometrically distributed somewhere 
within the dust of its host galaxy. For this example we 
choose a dust reddening law of /3dust = 2.8 (i.e. Ry — 
1.8). 

From these components we then simulate the final color 
and luminosities of the mock SNe la, and these are shown 
in Figure!?} Now we can see that a single monolithic color 
law fittedto this data will fit a line (solid black line) 
with slope intermediate between the input dust (shallow 
red dotted line) and SN intrinsic (steep blue dotted line) 
color laws. The luminosity correction error then depends 
on which side of the line the SN falls. Because dust con- 
tent is greater (on average) in more massive SN la host 
galaxies, more of these SNe la (magenta triangles) will 
fall above the line and will be over-corrected for dust 
(and vice versa for SNe la in low mass galaxies, green 
diamonds) . Similarly, if we split the sample in host mass 
and fit separate color laws, the high mass hosted SNe la 
have a lower /3, closer to the input dust law, while low 
mass hosted SNe la have a higher f3, closer to the input 
SN la intrinsic color law. T his result in consistent qual- 
itatively with that found by Sullivan et al. (2010) when 
splitting the SNLS SN sample by host mass. 

This set of input parameters (/?dust and /3sn) was pur- 
posely chosen for this example in order to produce the 
desired sign of the SN la host bias. A more detailed ex- 
ploration of this two-color effect across the whole space 
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Figure 7. Example Monte-Carlo simulation of competing SN la 
color variation of the SN la luminosity (steep dotted blue line) 
and foreground reddening by dust (shallow dotted red line). A 
fitted monolithic color law (black line) fits a slope between the two 
input slopes, while a color law fit split by host mass fits a color 
law (shallow dashed magenta line) for SNe la in high-mass hosts 
(filled magenta triangles) closer to the dust color law and a color 
law (steep dashed green line) for SNe la in low-mass hosts (open 
green diamonds) that is closer to the intrinsic SN la color law. All 
color law lines have been fixed to a common intercept to better 
illustrate the differing slopes. 



of input parameters will be presented in Paper III. For 
this work we focus on the trend of SN la Hubble residuals 
with host mass predicted to be produced by this effect. 
We thus generate a large sample (to reduce Monte Carlo 
noise) of mock SNe la using the above parameters, and 
the calculated Hubble residual bias curve is shown in 



Figure 11 of Section |5.5| (dotted magenta curve). 

Despite the favorable qualitative results from this sim- 
ulation, this two-color effect failed to produce sufficiently 
large offsets between low- and high-mass-hosted SN la 
Hubble residuals to match observations. Even the steep 
SN la color law used here produced a step size of about 
0.03 mag. Similarly, the struc ture does not match that 
observed in the data (Section |4.3[ ), particularly at high 
host masses where the dust-free elliptical hosts dominate 
and the dust over-correction effect disappears. However, 
while this effect alone cannot explain the SN la host bias, 
the correlation of dust content with host mass is certain 
to be impacting the observed colors of SNe la and may 
still be at least contributing to the observed host bias. 

5.2. Dependence of SN la Color on Progenitor 
Metallicity 

SNfactory data show a strong correlation of SN la light 
curve color (SALT2 c) with host galaxy gas-phase metal- 
licity, stronger than the analogous trend with host mass. 
This indicates a possible trend of SN la intrinsic color 
with metallicity which could easily drive some of the ob- 
served SN la Hubble residual bias. 

We quantify this color- metallicity trend here as follows. 
Because dust reddening in galaxies scales with mass and 
thus metallicity, we first isolate the subset of SNe la 
with extremely low reddening. We binned the data in 
bins of host metallicity of width 0.4 dex in the range 
7.7 < 12+log(0/7J) < 9.3, then used the aforementioned 
dust models (see Section 5.1) to calculate what fraction 



Figure 8. SALT2 color c vs. host galaxy metallicity for SNfactory 
SN la hosts. In bins of metallicity 0.4 dex wide, we fit for the 
weighted mean color (large green squares) of SN la hosts (small 
filled blue circles) after the rejection of extreme red outliers (small 
open red squares). The best fit trend fitted from these binned 
averages is shown as the thick magenta line. 



of SNe la in that bin would be reddened by more than 
0.05 mag. The lowest two metallicity bins had no outliers 
rejected (less than 5% were expected to have reddening 
above our cut) while the highest metallicity bin had 65% 
of its SNe la excised. In Figure [SJ we plot the SN la color 
(SALT2 c) versus host metallicity, with those hosts pass- 
ing the outlier rejection shown as small filled blue circles 
and those which failed this cut as open red squares. 

We note here that three SNe la in our sample were 
identified as showin g signatures of unburned carbon 
( Thomas et al.||2011| ), and these SNe are designated by 
large cyan circles (note none of these failed the outlier 



Thomas et al. (2011) noted that these SNe la 



cut). 

showed exceptionally blue colors, and indeed two of them 
appear significantly bluer than the color-metallicity lo- 
cus. Since they do not significantly pull the trend we 
wish to fit, they are left in our sample, but we note that 
this subclass of SNe la may not obey the apparent color- 
metallicity relation we have identified here in our low- 
extinction subsample. 

We use the sample of SNe la deemed to be "low- 
reddening" as described above to derive a linear fit to 
the trend of SN la color versus host metallicity, which is 
shown as the thick magenta line in Figure [8] The slope 
of the fitted trend is 0.028 ±0.010 mag/dex, such that an 
SN la with 1 dex higher metallicity would be 0.028 mag- 
nitudes redder (in SALT2 c, which scales directly with 
observed SN la B — V colors). 

This trend of SN la intrinsic color with metallicity may 
lend some credence to theoretical explorations of the ef- 
fect of metallicity on SN la color. Our measurement 
of 0.028 ± 0.010 mag/dex compares w ell to the value of 



0.025 mag/dex given by the models of Dominguez et al 



(2001|) over a co mparable range in m etallicity. Hoetiich 



et~al( ( |1998[ ) and |Sauer et aiTf ( |2008[ ) found that an in- 



crease m the progenitor metallicity for an SN la tended 
to yield redder colors in the o ptical, consis t ent w ith the 
trend we observe here, though Lentz et al. (20001 found 
a contradictory trend in their models. We note that our 
light curve colors were derived from fits to rest frame 
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BVR photometry, so our colors are not sensitive to the 
UV which may b e more strongly affected by metallicity 
than the opt ical dHoeflich et al.)|1998l |Lentz et al.||2000 



Sauer et al.||2008[ ). Interestingly, |Foley| ( 12012[ ) recently 



reported a correlation of the calcium velocity with host 
galaxy mass for a sample of SNe la. If their previously 
stated assumption that the calciu m velocity correlates 
directly with intrinsic SN la color ( Foley fc KasenpOll I 
holds true, then this would agree qualitatively with our 
color-metallicity trend measured from SNfactory data. 

Let us revisit the step in Hubble residuals found when 
splitting the sample into low and high metallicity hosts 
at 12 + log(0/H) = 8.8. The mean metallicity of the high 
metallicity hosts is (12 + log(0/H)} high = 9.02, and the 
mean of the low metallicity hosts is (12 + log(0/H)) ;ou; = 
8.35. Given the trend just fitted, this implies that the 
intrinsic colors of high-metallicity-hosted SNe la are on 
average 0.019 mag redder than those in low-metallicity- 
hosted SNe la (for the distribution of metallicities in our 
sample) . 

This possible intrinsic color discrepancy between 
SNe la in hosts of different metallicities may provide an 
explanation for the observed bias of corrected SN la lu- 
minosities with host mass (and metallicity). Let us con- 
sider the simplest scenario where the stretch-corrected 
peak luminosity of SNe la is independent of metallicity, 
but their intrinsic colors depend on metallicity. SN la 
cosmology analyses typically use a color-luminosity cor- 
rection factor of about j3 ~ 3, which implies that the 
luminosities of SNe la in high-metallicity hosts may be 
over-corrected for their color by about Am = 0.057 mag. 
The fact that the observed Hubble residual step between 
high- and low-metallicity hosted SNe la is larger than this 
could indicate that other effects may be contributing to 
the observed SN la host bias. 

The full picture of SN la intrinsic color and its implica- 
tions for color correction techniques is likely to be more 
complic ated , including such effects as were outlined in 
Section 15.11 Our main interest is the structure of SN la 
Hubble residuals as a function of host ma ss produced by 
this trend, which we show in Figure 11 (dashed green 
curve). This model does not accurately trace the full 
structure of the observed host bias, but interestingly this 
model predicts a flattening of the Hubble residual trend 
at high host mass (log(M*/M ) ~ 11.0) similar to that 
seen in the data. This is due to the flattening of the 
galaxy mass-metallicity relation at this mass scale, and 
may point to the importance of metallicity in the struc- 
ture of the SN la Hubble residual host trend. 

5.3. Young versus Old SN la Progenitor Ages 

While the interplay of SN la intrinsic colors and dust 
are intriguing candidates for explaining the observed host 
bias, another potential source of the bias may arise from 
SN la progenitor ages. Studies have shown that the SN la 
rate in a given galaxy depends on both its stellar mass 
and star- format ion rate, and we show here that the way 
those quantities are distributed along the galaxy mass se- 
quence has direct implications for SN la progenitor ages. 

The distribution of stellar mass along the g alaxy mass 
sequen ce is well-described by a (modified) |Schechter 
( 1976 ) function, and was first compared to the SN la 
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Figure 9. Top: Histogram (solid blue line) and generalized his- 
togram (dashed blue line) of SN la host galaxy masses from the 
full sample of SNe la discovered by the SNfactory search, compared 
to a fit (dashed red curve) composed of a component proportional 
the distribut ion of stellar mass (dashed green curve) in the lo- 
cal universe ( |Schechter||1976| |BeIl et al ||2003| [B 
and a component proportional to the distributii 



IBorch et al. 



[20CJ6] | 



tion of star torma- 



componcnt proportional 
tion (dashed magenta curve) in the local universe l |Salim et al.| 
2007 . El baz et al.|2007) . Bottom: Ratio of young SN la progeni- 
tors (magenta curve in top panel) to all SNe la (red curve in top 
panel) as a function of mass. SNe la undergo a transition to older 
progenitors at about 1O 1O ' 5 M0. 



massive galaxies are more int ensely star-forming than 



their high-m ass counterparts (Salim et al. 2007 



Elbaz 



et al. 20071, making the distribution of star-formation 



along the galaxy mass sequence peak at lower masses 
and have a stronger low-mass tail. Following the SN la 
rate dependence on these two quantities, the SN la host 
mass distributio n can be fit as a linear combinatio n of 
the stellar mass (Bell et al. 2003;|Borch et al.||2006|) and 



star- formation ([Salim ct al. 20071 |tilbaz et al.puuf ) dis- 
tributions as measured from SUSS^ 

We show in the top panel of Figure [9] the results of 
fitting the host mass distribution for the entire sample 
of SNe la discovered by SNfactory with these two com- 
ponents. We note here that the only free parameter in 
this fit is the relative ratio of the two components, as 
their shapes and positions along the galaxy mass scale 
are fixed by measurements of large external galaxy sam- 
ples. This fit yielded a measurement of the fraction of 
SNe la arising from young progenitor systems (i.e. the 
"prompt" fraction) of 67 ±6%. This value was consistent 
with that derived from SN la rate coefficients from other 



studies (Scannapieco & Bildsten |2005j ?). 



host luminosity function in Gallagher et al. (2005). Less 



In the lower panel of the same Figure we show the 
prompt fraction as a function of host galaxy mass, which 
shows that the fraction of SNe la from young progen- 
itors remains very high (> 90%) until it undergoes a 
smooth transition to being subdominant over the mass 
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range 10.0 < log(M*/M ) < 11.0. This shows that 
SNe la in low mass hosts originate primarily from young 
("prompt") progenitors while the majority of SNe la 
in high mass hosts come from older ("tardy") progen- 
itors, with the equality point being reached at about 
log(M*/M0) = 10.5. We note that this coarse two- 
component model would be better investigated with a 
full SN la delay time distribution and its implication for 
SN la ages along the galaxy mass sequence, and will be 
presented in the Paper IV. 

Thus we see that SNe la undergo a significant change 
in progenitor age at a galaxy mass scale that is close to 
scale at which corrected SN la Hubble residuals undergo 
a transition. This could be a potential explanation of 
the observed host bias in SN la Hubble residuals if the 
luminosity of SNe la from young progenitors is different 
from that of old progenitors, or alternatively if the mean 
colors of these two subsets of SNe la differ. This could 
be a dependence of SN la luminosity (or color) on age, 
or it could be indicative of two different SN la progenitor 
channels which operate on different timescales and have 
different mean luminosities. 

If the prompt/tardy transition is the driver of the host 
bias, we might expect the shape of the prompt frac- 
tion transition to reflect that seen in the Hubble residual 
data. We thus fit the prompt fraction model to the ob- 
served Hubble residual data by fitting for a mean prompt 
Hubble residual value and a mean tardy Hubble residual 
value. Minimizing the x 2 °f this model yields a best 
fit luminosity difference between the prompt and tardy 
SNe la of 0.15 mag. This best fit model is shown in Fig- 
ure |11| The steepness of the prompt fraction transition 
appears slightly shallower than the observed structure in 
the SN la Hubble residual trend with host mass, so this 
model may be too simple. However, this trend is closer 
to the observed trend than the other models we consid- 
ered, so we examine some important implications if it is 
true. 

This transition host mass for SN la progenitor ages 
is perhaps unsurprising given that the normal galaxy 
population undergoes a transition from primarily ac- 
tive to pass ive galaxies at th is sa me mass scale . From 



SDSS data, |Bell et al.| fl2003| ) and |Baldry et aL| ( [2004 1 
fit the s tellar mas s func tions - which are well fit by 
classical Schechter (1976) functions - for red (passive) 
and blue (active) galaxies. Both studies found that 
the number density (and thus mass density) of red and 
blue galaxies become equal at a galaxy mass scale of 
about log(M*/M ) = 10.3. A similar transition between 
red and blue galax i es is also evident at higher redshifts 



(iBorch et al.l|2006| |Bundy et al.||2006| |Brammer et al 
2011 1 , with the transition mass increasing as redshift in- 



creases. 

If the SN la host bias is indeed driven by progen- 
itor age, then the evolution in redshift of this galaxy 
color (and thus age) transition mass has important con- 
sequences for SN la luminosity corrections. For exam- 
ple, if SN la host masses are used to correct luminosities 
via a stepwise correction where luminosities of SNe la 
in hosts below some specific mass are adjusted with a 
constant shift, then the evolution of the galaxy age tran- 
sition mass might imply that the SN la threshold host 
mass must necessarily evolve in redshift to avoid bias. 

If progenitor age is indeed the cause of the observed 



host bias, then failure to correct SN la luminosities will 
result in cosmological biases due to the evolution of 
mean SN la progenitor ages. Suppose for simplicity that 
the luminosity offset between high- and low-mass hosted 
SNe la implies that young (prompt) SNe la are intrinsi- 
cally fainter than old (tardy) SNe la, then the shift in the 
relative fraction of these two subsamples between high 
and low redshift will mean that the average corrected 
luminosities of SNe la differ at high and low redshift. 
Using th e mean mass-SFR relation at high-redshift mea- 
sured by Elbaz et al. (2007) adjusted to the same IMF as 
ours, we calculate that by z ~ 1 the fraction of prompt 
SNe la should increase to 95% (this is similar to the re- 
2007] 



suit of Howell et al 

of low stretch, i.e. "tardy 



who noted that the fraction 
SNe la should decrease from 
about 25% to about 1% beyond redshift z ~ 1). Thus 
if the two sub-populations of SNe la differ in luminosity 
by 0.15, then an increase of 25% of the prompt fraction 
implies SNe la at high redshift will be about 0.038 mag 
fainter than those at low redshift after stretch and color 
corrections. 

5.4. Metallicity Dependenc e of SN la 56 Ni Yie ld and 



Rise Time from Bravo et al. (2010) 



The correlation of galaxy mass with metallicity has 
made metallicity a leading contender for explaining the 
bias of SN la Hubble residuals with host mass. In SN la 
models, metal licity affects the opacity and thus the re - 
sulting colors (Domfnguez et al. 2001; Kascn et al. 2009), 
and th e dynam ics leading to explosion may also be af- 



fected (Townsley et al. 2009). The largest effect may 
arise from the predicted decrease in the production of 
56 Ni due to an increa se in neutronizatio n with meta llicity 



( |Timmes et~aT] ; 2003). In related work, |Piro| ( |2008[) have 
suggested that additional neutronization occurs during 
the "simmering" phase, such that metallicity-induced 
neutronization will not become apparent until metallici- 
ties above 2/3Zq axe reached. 

Rather than attempt to transfo rm our data into mo del 
space, we employ the results of Bravo et al.| ( |2010[ ) to 
determine relative brightnesses after stretch correction 
as a functi on of metallicity pr edicted by theoretical cal- 
culations. |Bravo et al. (2010) provide formulae for pre- 
dicted peak brightness and light curve width, allowing a 
stretch correction to be applied to th e models, as with 
the observations. |Bravo et al. (2010) account for burn- 
ing that does not reach nuclear statistical equilibrium 
and find that neutronization during the simmering phase 
does not strongly affect the relative brightness trend with 
metallicity. They calculate models with and without a 
metallicity dependence for the deflagration to detonation 
transition density pddt] the former model produces a 
stronger de pend ence on met allicity. 

In Figure 10 we plot both |Bravo et aL ( 20101) models 
(note the expanded y-axis scale). While the JBravo et al. 
(20101 models make predictions for bolometric quanti 



ties, we have found with SNfactory data (Scalzo et al., 
2013, in prep.) a linear mapping between the B-band 
and bolometric quantities. This implies that the .B-band 
and bolometric model predictio ns should be simila r and 
it is worthwhile to compare the |Bravo et al. (2010) pre- 
dictions directly with the -B-band Hubble residuals. It 
is quickly apparent that both models predict a trend of 
corrected luminosity with host mass that is opposite to 
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Figure 10. Binned avera ge Hubble re siduals, as in Figure|6] com- 
pared to the two models of Bravo ct al. (2010): the constant Pddt 
model (dashed-dotted cyan curve] and the metallicity dependent 
Pddt model (thin solid cyan curve). For reference we also show 
the Hubble residual averages for the two "plateau" regions as the 
thick red lines. 



that observed in the data. This is due to the dimming of 
SNe la from the neutronization effect at high metallicity 
having a larger effect on the final SN la luminosity than 
stretch. 



Interestingly, Kasen et al. (2009) predi ct an effect 



whose sign is opposite to that predicted by |Bravo et al. 



(2010) and coincides with the sign of bias seen in the 
data, as was previously noted by Sullivan et al. ( 2010 ). In 
those more detailed models the metallicity has a stronger 
effect on the lightcurvc stretch than on the luminosity via 
the effect of temperature on opacity. Due to the compu- 
tational demands of these models the sampling in metal- 
licity is too coarse to make a quantitative comparison 
with our data at this time. 

5.5. Summary and Comparison of Models 

We summarize the potential drivers of SN la host bias 
presented here. In FigurefTTlwe show the predicted shape 
of the mean SN la Hubble residuals versus host galaxy 
mass for the models presented above, as well as the best 
fit linear trend from Section [4] 

For the sake of comparison, we show a model repre- 
sented by the two plateau values (i.e. constant Hub- 
ble residual value) outside the transition mass region of 
9.8 < log(M*/M Q ) < 10.4 and connected with an error 
function. This functional form is exactly what would be 
expected for a step function convolved with a Gaussian, 
and our best fit for this functional form yielded a step 
function which transitions at log(M*/JW0) = 10.21 con- 
volved with a Gaussian of width a = 0.47 dex. This 
model, shown in Figure [TT] as the thin solid red line, pro- 
vides the best fit to the data. We note that the Gaussian 
width of ct = 0.47 dex for this model is much larger 
than the typical measurement error of 0.1-0.2 dex for 
host masses, and indeed exceeds even the most conser- 
vative estimates of systematic stellar mass uncertainties 
(see Section pi). The driver of this transition width may 
indicate an astrophysical effect such as the width of the 
galaxy mass-metallicity relation (in mass) or the charac- 
teristic width of the blue-red galaxy transition (e.g. Sec- 
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tion 5.3). Greater statistics of SN la host galaxy metal- 
licities or ages will likely be necessary to show whether 
the Hubble residual transition is sharper with respect to 
either of these properties. 

In Table |3| we present the % 2 of the six models plotted 
in Figure [TTJ when compared to the full sample of 601 
SNe la from the four data sets, as well as the improve- 
ment in x 2 f° r each model and the number of degrees 
of freedom n^of ■ Note that five cosmological parameters 
were applied to the data so that reduces the starting n^of 
to 596. For reference, we also report the x 2 if no trend 
is considered, which is exactly the \ 2 °f the standard 
cosmological fit for these SNe. 

Our first conclusion is that the 2-color and color- 
metallicity models are unlikely to be the sole driver of 
the observed host bias. However, the trend of SN la red- 
dening by foreground host galaxy dust is based on obser- 
vationally determin ed galaxy properties ( Lee et al.|2009 



Garn fc BestpOlO ), and SN la colors are probably being 
affected by this trend to some degree. If so, the magni- 
tude of the effect depends on the ability to successfully 
recover the dust reddening for each SN, which in turns 
depends on knowledge of the intrinsic colors of SNe la 
and the true dust extinction law. The color-metallicity 
trend agrees well wit h the oretical models and some ob- 
servations (see Section 5.2 ) so may still be affecting SN la 
color corrections. 

Our second conclusion from this analysis is that the 
SN la progenitor age model presents the best quanti- 
tative agreement with the observed host bias of all the 
physical models. This underscores the importance of 
considering age as a cont ender for driving th e host bias. 
Thirdly we found that the Bravo et al. ( 2010 ) neutroniza- 
tion model exhibits poor agreement with the SN la host 
bias structure, as it predicts a trend opposite to that 
observed in the data. Finally we note that the analyti- 
cal error function model yielded the best fit to the data. 
Along with the favorable agreement of the prompt /tardy 
model, this result may indicate that the SN la Hubble 
residuals are better represented by a bimodality with re- 
spect to their host properties rather than a continuous 
linear trend. 

We note that other scenarios may exist that could pro- 
vide an explanation for the observed SN la host bias. For 
example, if the absorption properties of dust (i.e. Ry) 
vary with metallicity (a troubling possibility for stan- 
dard candle cosmology), this could provide a less ortho- 
dox explanation of the observed trend of SN la (3 with 
host mass, and would ce rtainly result in bia sed SN la 
color corrections (but see Chotard et al.|[20lT ). We have 
thus far neglected possible effects originating from red- 
dening of SNe la by some circum-stellar medium (CSM) 
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Figure 11. Binned average Hubble residuals, as in Figure^] compared to predictions for several models: (i) simple two-color model of 
Section pTl] (dotted magenta curve); (ii) trend of SN la intrinsic color with metallicity as discussed in Section |5 . 2| (dashed green curve); (iii) 
model where SNe la originate fro m "prompt" and "tardy" populations whose luminosities differ by the observed Hubble residual plateau 
offset, as described in Section |5,3| and the "prompt fraction" of SNe la changes along the galaxy mass sequence; (iv) a simple linear trend 
fitted to the data (thick orange clashed line); and (v) low-mass and high-mass Hubble residual plateaus connected by an error function, 
which is exactly a step function convolved with a Gaussian (thin solid red lines). 



deposited in the immediate SN vicinity by material shed 
in the later st ages of the progenitor system evolution 
( |Goobar|p008l ). If the CSM has a different R v than 
the ISM dust, especially if CSM dust content depends 
on progenitor metallicity, then this could produce simi- 
lar color correction er rors as those predicted in our two 
color model of Section [O] These and other explanations 
are worth investigation, but detailed discussions of their 
implications are beyond the scope of this work. 

In summary, we demonstrated that dust, age, and 
metallicity all vary along the galaxy mass sequence. 
While dust effects alone could not quantitatively explain 
the observed SN la host bias, dust is probably still affect- 
ing SN la luminosity corrections in a small but important 
way. Metallicity-based effects show a characteristic flat- 
tening of Hubble residuals at high galaxy masses, due 
to the flattening of the galaxy mass-metallicity relation, 
which matches the observed SN la host bias structure. 
This may point to metallicity as a leading contender for 
driving the host bias. However we also showed that the 
likely evolution of SN la progenitor age along the galaxy 
mass sequence also shows favorable agreement with the 
data, so age should not be ruled out as a possible driver 
of the host bias. 

None of these effects acting alone seem able to mimic a 



transition as steep as that observed in the binned SN la 
Hubble residuals. It may be the case that the age, metal- 
licity, and dust trends are all present in some combi- 
nation. We believe that discrimination between these 
properties will require either (or probably both) a more 
detailed study of the immediate SN la environments or 
a complete accounting of intrinsic SN la color variation 
via constraints on host galaxy dust reddening. A larger 
sample of measured SN la host ages and metallicities may 
elucidate whether the Hubble residual host bias is more 
strongly correlated with these properties than with host 
mass. 

6. CORRECTING HOST BIAS IN SN la DATA 

In light of the observed SN la luminosity host bias, 
a major concern for future SN la cosmological stud- 
ies is how to effectively correct for this effect in SN la 
data. Currently, the observed correlation of SN la Hub- 
ble residuals with their host galaxy masses has motivated 
the possible use of SN la host mass as a third brightness 
correction parameter (in addition to stretch and color) 
in SN la cosmological analyses. Indeed, some authors 
have already performed SN la cosmo logy fits using a lin- 
ear host-based brightness correction ( Sullivan et al.|20TT I 
or a step-wise brightness correction for host mass where 
brightnesses of low-mass-hosted SNe la are adjusted by 
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a constant amount ( Suzuki et al 



2012). 



In this Section we briefly describe the potential contri- 
bution of host corrections to the SN la cosmology error 
budget, as well as potential sources of bias in SN la cos- 
mology if host corrections are not used. In particular, 
we seek to estimate two quantities: (i) the amount of er- 
ror introduced to SN la Hubble residuals from the host 
correction itself, and (ii) the bias in mean luminosity of 
high-redshift (z ~ 1) SNe la if host corrections are not 
applied. 

Metallicity is currently perceived to be a leading can- 
didate for explaining the origin of the SN la host bias, 
but we showed above that dust and progenitor age may 
be driving part or all of the observed host bias. We will 
focus on dust effects in Paper III, and age effects in Pa- 
per IV, but here we explore the potential implications if 
the trend is metallicity-driven. 

Specifically, we consider a model where metallicity af- 
fects the luminosity of SNe la in a linear way with a slope 
equal to that measured from SNfactory data in Section 2j 
and assume that SN la luminosity corrections will be de^ 
rived from photometrically-estimated host galaxy stellar 
masses. The major sources of systematic errors under 
this scheme would be the systematics in galaxy stellar 
population synthesis (SPS) techniques to estimate host 
masses, the dispersion in (and evolution of) the galaxy 
MZ relation, and the difference between the inferred aver- 
age metallicity of the SN la host galaxy and the metallic- 
ity of the SN progenitor itself. We address and quantify 
each of these effects in turn. 

6.1. SPS systematics 

If host galaxy properties are used to correct SN la lu- 
minosities, then the SN la error budget will be increased 
by the systematics involved in this process. In partic- 
ular, if host masses are used as correction factors, then 
the systematics accompanying the application of stellar 
population synthesis (SPS) techniques will impact SN la 
cosmology. Though a census of the full body of SPS sys- 
tematics is beyond the scope of t his work (see thor ough 
discu ssions in, e.g., Gallazzi & Bell 2009; |Conroy et aL[ 
|2009| regarding mass-to- light ratio recovery efficiencies 
and the impact of stellar models on fitted mass-to-light 
ratios), we summarize here a few key elements that will 
s urely impact host-bas ed SN la luminosity corrections. 



Conroy et al. (2009) showed that uncertainties in the 
contribution of 1P-AGB stars to galaxy luminosities and 
uncertainties in the shape of the stellar initial mass func- 
tion (IMF) can result in stellar mass uncertainties of 
order 0.3 dex. The ability to recover mass-to-light ra- 
tios is limited by how well photometry can constrain 
the galaxy star formation history (SFH). We showed in 
Paper I that full photometric coverage from UV to the 
NIR can provide mass-to-light ratios with a dispersion of 
about 0.1 dex, but even a single color constrains these to 
about 0.3 dex. Gallazzi & Bell (2009) showed that erro- 
neous assumptions about priors on galaxy SFHs, such as 
the frequency and intensity of bursts of star formation, 
can bias mass-to-light ratio values by as much as about 
0.2 dex. 

Considering these effects together, galaxy SPS tech- 
niques carry with them a systematics floor that may be 
as high as 0.4 dex in stellar mass. We note that this 
is a conservatively high value for SPS uncertainties if 



one assumes that all galaxy masses are affected by TP- 
AGB, SFH, and IMF uncertainties at the maximum level 
estimated by previous authors. If the stellar mass un- 
certainty were truly this high, then it would encompass 
almost the entirety of the 0.47 dex width of th e step func- 
tion Hubble residual model from Section [53] This would 
require a physical model of SN la Hubble residuals that is 
a true step function in host mass, which is very unlikely. 
Typical host mass uncertainties for the SNfactory sample 
are at the 0.15 dex level and are driven by SFH uncer- 
tainty (Section [51 also Paper I), but we calculate this 
maximal host mass error of 0.40 dex to demonstrate the 
maximal cosmological uncertainty host mass estimates 
might introduce. 

Propagating this uncertainty in the host mass esti- 
mate through our fitted trend of SN la luminosities with 
mass implies that the resultant uncertainty in the host- 
correction terms for SN la luminosities will have an ad- 
ditional scatter of about 0.013 mag. We note than any 
coherent error in the estimation of host galaxy stellar 
masses would be absorbed into the standardized SN la 
luminosity (i.e. Mb) so long as it is applied consistently 
across all redshifts, and thus would not bias the measure- 
ment of cosmological parameters. 

6.2. Dispersion in (and evolution of) the MZ relation 

The use of photometric SN la host galaxy masses as 
proxies to correct for a metallicity-driven effect will have 
some systematic errors resulting from the dispersion in 
the galaxy MZ relation. At high metallicities, this disper- 
sion is as low as 0.1 dex, but increases to about 0.3 dex 



at low metallicities ( Tremonti et al. [2004 ) . We use the 
observed trend of SN la Hubble residuals with metallic- 
ity from SNfactory (Section [3]) to calculate that, even 
at the extreme high MZ dispersion of 0.3 dex, the MZ 
relation dispersion would introduce a dispersion in host- 
mass-corrected SN la luminosities of about 0.032 mag. 

If indeed the observed host bias is driven by progen- 
itor metallicity, then the use of SN la luminosity cor- 
rections derived from host galaxy masses must account 
for evolution of the MZ relation with redshift in order 
to av oid biasing m e asure ments of cosmological parame- 
ters, 
lation 



Zahid et al. (2011) measured the galaxy MZ re- 
at z '■-' 0.8 from the DEEP2 survey and found 



that galaxies at high redshift have metallicities approxi 
mately 0.15 dex lower than their counterparts in the local 
universe. Again using the host bias trends measured in 
Section [3[ this would imply a systematic bias of about 
0.016 mag between the corrected luminosities of low- and 
high-redshift SNe la. Though very small, this bias would 
be present if SN la host masses are used to correct a 
metallicity-driven SN la luminosity effect without con- 
verting to the redshift-dependent metallicity implied by 
the evolving MZ relation. 

Finally we consider the potential bias in SN la lumi- 
nosities if host-based corrections are not applied. The 
main driver of this bias is the difference in average host 
galaxy mass between the low- and high-redshift SN la 
samples. For example, the mean host mass and redshift 
of the low redshift SN la sample from Kelly et al. (2010) 
0.035 and log (M«/ M M ) = 10.67, wEIeT oFTne 



are z 



high-redshift SNLS sample of |Sullivan et al.| ( |2010[ ) these 
are z = 0.634 and log(M*/iW ) = 9.79. UsmgHie MZ 
relation of T04 and a MZ offset of 0.1 dex (scaled down 
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approximately by redshift) due to its evolution with red- 
shift, this implies that current samples of high-redshift 
SN la hosts are likely to be 0.29 dex more metal poor 
than the low-redshift sample. This would imply an off- 
set in stretch- and color-corrected SN la luminosities of 
0.031 mag in current SN la samples. However for SNfac- 
tory data at low redshift (with median stellar mass of 
log(M st /Af Q ) = 10.1), this implies the offset between 
high and low redshift SN la host metallicities would be 
about 0.2 dex, corresponding to a magnitude offset of 
0.021 mag. 

6.3. Discrepancies between host and progenitor 
metallicities 

We now turn to the systematic error on host-corrected 
SN la luminosities introduced by the difference between 
the measured (or inferred) host metallicity and the 
metallicity of the SN la progenitor. An individual galaxy 
has an internal distribution of stellar metallicities with 
non-negligible dispersion, often with gradients as high 
as 0.5 dex per scale radius (e.g. Spolaor et al. 20101, 
and the SN progenitor metallicity is drawn at random 
from this distribution. Based on the metallicity trend 
observed above with SNfactory data (Sectionp?]), a galaxy 
with internal metallicity dispersion of 0.5 dex would in- 
troduce a random error on SN la brightnesses of about 
0.053 mag, which is much smaller than the observed dis- 
persion in corrected SN la luminosities (~ 0.15 mag) or 
the offset between the corrected luminosities of high- and 
low-metallicity- hosted SNe la (~ 0.1 mag). Thus despite 
the potential differences between measured SN la host 
metallicities and the progenitor metallicities themselves, 
these would ultimately have only small effects on the 
corrected brightnesses of SNe la. The range of observed 
galaxy metallicities has a much larger span than the typ- 
ical internal metallicity dispersion of a single galaxy, so 
ultimately such a correction would introduce an error 
whose magnitude is smaller than the bias it is intended 
to correct. 

It is also important to consider potential biases intro- 
duced by the difference between the SN la progenitor 
metallicity at the time of its formation and the metallic- 
ity of the host galaxy measured at the time of explosion. 
Chemical evolution in the host galaxy of a SN la means 
that the average current metallicity may be significantly 
higher than it was at the time the SN la progenitor sys- 
tem was formed (especially for long del ay times of several 
Gyr) , by as much as perhaps 0.6 dex ( |Bravo fc Badenes 
20111. If all SNe la have the same (average) offset be 



Table 4 

SN la Host Correction Errors 



tween progenitor and host metallicity, then this offset 
could be absorbed into the standardized luminosity Mb ■ 
More specifically, so long as the average offset between 
host and progenitor metallicity remains constant at all 
redshifts, then the mean Mb would be valid for all red- 
shifts. 

However, an evolution of the average SN-host metal- 
licity offset would produce a bias between high- and low- 
redshift samples. This would most likely arise if the av- 
erage delay time of SNe la evolves with redshift, which 
we expect it should. To explore this in a simple fash- 
ion, let us consider the simple segregation of SNe la into 
short delay times ( "prompt" SNe la) and long delay times 
( "tardy" SNe la) and assume that all prompt SN la pro- 
genitors are at the same metallicity as their host whereas 
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- SFH Uncertainty 
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0.007 
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0.011 
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0.013 



MZ Relation 

- Mass-based Z-correction 

- MZ Evolution (vs. z) 
SN-Host Z Mismatch 

- Host Z dispersion 

- Delay time Z evolution 



0.016 



0.016 



0.032 



0.053 



all tardy SN la progenitors are actually 0.6 dex lower 
metallicity than their hosts. In this scheme the average 
SN-host metallicity offset will change from low- to high- 
redshift proportionally to the fra ction of tardy SNe la. 
As discussed above in Section [53} we calculated in Chil- 
dress et al. (2012c, in prep.) that the current prompt 
fraction at low redshift is about 70%, while at high red- 
shift it should increase to 95% (due to the increasing 
intensity of star-formation) . Thus the tardy fraction de- 
creases by 25%, so the average offset between SN and 
host metallicity would decrease by 0.15 dex. This means 
the average SN la luminosity zeropoint at high redshift 
would be brighter by about 0.016 mag. 

6.4. Summary and Outlook 

In this Section we investigated possible contributions 
of host galaxy systematics to the scatter in SN la Hub- 
ble residuals if hosts are used as a third luminosity cor- 
rection parameter, as well as the bias in high-redshift 
SN la luminosities if a host correction is not applied. We 
proceeded under the simple assumption that metallicity 
linearly affects the SN la luminosity, and calculated bias 
and dispersion values which are summarized in Table 111 

The sources of dispersion in corrected SN la luminosi- 
ties generally arise from systematic uncertainties in the 
estimation of host galaxy properties and the differences 
between global host metallicity and the metallicity of the 
SN progenitor itself. Even in the most pessimistic sce- 
nario, these effects have a quadrature sum of 0.063 mag, 
which is significantly smaller than the observed scatter 
in SN la Hubble residuals. Conversely, the failure to ap- 
ply host corrections to SN la luminosities will bias SN la 
Hubble residuals at high redshift due to the evolution of 
mean progenitor metallicity at the level of 0.031 mag. 
We thus claim that host corrections should be applied 
for current SN la samples as these will correct a critical 
source of bias while contributing negligibly to the SN la 
Hubble residual scatter. 

Finally, we caution that a proper accounting of the evo- 
lutionary effects of galaxy metallicity with redshift is crit- 
ical to avoiding bias in SN la host corrections. Failure to 
incorporate a redshift-dependent galaxy mass-metallicity 
relation would bias high redshift SN la Hubble residuals 
by 0.016 mag. A similar amount of bias could be intro- 
duced if the mean delay time between the SN la pro- 
genitor production and its explosion causes a systematic 
offset between host and SN metallicity which is not ac- 
counted for. The former effect is well studied and can be 
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directly compensated for, while the latter is dependent 
on the SN la delay time distribution which is currently 
still a topic of vigorous study. 

7. CONCLUSIONS 

Using SN la light curves and host galaxy properties 
from the Nearby Supernova Factory (SNfactory), we 
showed that the stretch- and color-corrected SN la lumi- 
nosities exhibit a bias with respect to the masses, specific 
star- formation rates, and metallicities of their host galax- 
ies. We also showed analogous trends exist for model-free 
observational quantities such as host absolute magnitude 
and color. These results are i n good agreement with 
the findings of previous authors (Kelly et al.|2010 |Lam 



illy ( 
peitl et al.i20T0||Suliivan et al.|20lO||konishi"et ai".J2Ull 



Gup ta et al.||2011[[D'Andrea et al.||201ip 
When we combine our results with 



those of 



et al fl2010| ), |Sullivan et aL| ( |2010[ ), and |Oupta et ah 
(|2U11 ), we demonstrate that stretch- and color-corrected 
SIN la Hubble residuals of SNe la in high- and low- 
mass hosts differ by 0.077 ± 0.014 mag, a 5.6cr detec- 
tion of SN la host bias. Furthermore, we show that 
when Hubble residuals for these data sets are derived 
with common SN la luminosity correction parameters 
(a, P) the shape the Hubble residuals trend appears to 
plateau at low and high host masses, with a difference 
of 0.086 ± 0.016 mag that transitions quickly in the host 
mass range 9.8 < log(M*/M ) < 10.4 (see Section [O 



and Figure [6]). 

We then examined several physical effects which could 
potentially drive the observed trend of SN la Hubble 
residuals with host galaxy properties. Most importantly 
we highlighted the fact that dust, age, and metallicity 
all vary along the galaxy mass sequence and could con- 
tribute to the SN la host bias. Dust alone could not ex- 
plain the entire bias observed, but both metallicity and 
age could produce Hubble residual trends with similar 
magnitude and structure to those observed in the data. 

Finally we considered the cosmological implications of 
using (or not using) SN la host properties to correct their 
luminosities. We showed that systematic errors in the 
estimation of SN la host galaxy properties would yield 
a maximum total dispersion of 0.06 mag, but failure to 
correct the observed Hubble residual host trends would 
bias the luminosities of high-redshift (z ~ 1) SNe la by 
about 0.02-0.04 mag. We thus concluded that host cor- 
rections in SN la cosmology contribute negligibly to the 
SN la error budget but correct a critical source of bias 
and are recommended to be used at this time. However 
we cautioned that a full understanding of the astrophys- 
ical origin of the bias is necessary in order to derive a 
correct redshift dependence of the host corrections. 

Correction of the observed SN la Hubble residual host 
bias is a critical step for the future of SN la cosmology. 
The obvious current means to account for this bias is 
to apply a third empirical SN la luminosity correction 
derived from host galaxy properties, but the ideal future 
goal would be to uncover the physical origin of this trend 
and find a more direct means of correcting it. We have 
shown here that many aspects of SN la progenitors and 
environments change along the galaxy mass scale and 
may be altering the observed luminosities and colors of 
SNe la through several mechanisms. We suggest that the 
preferred next step for investigating the SN la host bias 



is to build a sample of SNe la whose environmental prop- 
erties (particularly age, metallicity, and dust reddening) 
can be constrained on a local scale. The ideal hunting 
ground for such SNe la is in the low redshift universe 
where measurement of spatially resolved galaxy proper- 
ties is possible. 

Facilities: UH:2.2m, GALEX, Shane (Kast Double 
spectrograph), Keck:I (LRIS), Blanco, SOAR, Gem- 
ini:South 

Acknowledgments: Based in part on observations made with the 
NASA Galaxy Evolution Explorer. GALEX is operated fro NASA by 
the California Institute of Technology under NASA contract NAS5- 
98034. The authors graciously acknowledge support from GALEX 
Archival Research Grant #08-GALEX508-0008 for program GI5-047 
(PI: Aldering). This work was supported by the Director, Office of Sci- 
ence, Office of High Energy Physics, of the U.S. Department of Energy 
under Contract No. DE-AC02-05CH11231; the U.S. Department of 
Energy Scientific Discovery through Advanced Computing (SciDAC) 
program under Contract No. DE-FG02-06ER06-04; by a grant from 
the Gordon &: Betty Moore Foundation; in France by support from 
CNRS/IN2P3, CNRS/INSU, and PNC; in France by support from 
CNRS/IN2P3, CNRS/INSU, PNC, and Lyon Institute of Origins under 
grant ANR-10-LABX-66; and in Germany by the DFG through TRR33 
"The Dark Universe" . This research used resources of the National En- 
ergy Research Scientific Computing Center, which is supported by the 
Director, Office of Science, Office of Advanced Scientific Computing 
Research, of the U.S. Department of Energy under Contract No. DE- 
AC02-05CH11231. We thank them for a generous allocation of storage 
and computing time. HPWREN is funded by National Science Founda- 
tion Grant Number ANI-0087344, and the University of California, San 
Diego. The Centre for All-sky Astrophysics is an Australian Research 
Council Centre of Excellence, funded by grant CE1 10001020. 

The authors would like to thank the excellent technical and scientific 
staff at the many observatories where data was taken for this paper: 
the University of Hawaii 2.2m telescope, Lick Observatory, Keck Ob- 
servatory, the Blanco 4m telescope, the SOAR telescope, and Gemini 
South. Some data presented herein were obtained at the W. M. Keck 
Observatory, which is operated as a scientific partnership among the 
California Institute of Technology, the University of California, and the 
National Aeronautics and Space Administration; the Observatory was 
made possible by the generous financial support of the W. M. Keck 
Foundation. We wish to recognize and acknowledge the very signifi- 
cant cultural role and reverence that the summit of Mauna Kea has 
always had within the indigenous Hawaiian community, and we are ex- 
tremely grateful for the opportunity to conduct observations from this 
mountain. We also thank Dan Birchall for assistance with SNIFS ob- 
servations. We are very grateful to David Rubin for providing SALT2.2 
light curve fits to the CfA light curves in advance of the forthcoming 
Union3 analysis. We also thank Josh Meyers and Dan Kascn for use- 
ful discussions. We thank the anonymous referee who provided very 
helpful comments. 

Some of the data analyzed here were obtained from the Sloan Dig- 
ital Sky Survey Eight Data Release (SDSS-III DR8). Funding for 
SDSS-III has been provided by the Alfred P. Sloan Foundation, the 
Participating Institutions, the National Science Foundation, and the 
U.S. Department of Energy Office of Science. The SDSS-III web site 
is http://www.sdss3.org/ SDSS-III is managed by the Astrophysical 
Research Consortium for the Participating Institutions of the SDSS- 
III Collaboration including the University of Arizona, the Brazilian 
Participation Group, Brookhavcn National Laboratory, University of 
Cambridge, Carnegie Mellon University, University of Florida, the 
French Participation Group, the German Participation Group, Har- 
vard University, the Instituto de Astrofisica dc Canarias, the Michigan 
State/Notre Damc/JINA Participation Group, Johns Hopkins Univer- 
sity, Lawrence Berkeley National Laboratory, Max Planck Institute 
for Astrophysics, New Mexico State University, New York University, 
Ohio State University, Pennsylvania State University, University of 
Portsmouth, Princeton University, the Spanish Participation Group, 
University of Tokyo, University of Utah, Vandcrbilt University, Uni- 
versity of Virginia, University of Washington, and Yale University. 

REFERENCES 



Aihara, H., et al. 2011, ApJS, 193, 29 

Aldering, G., et al. 2002, in Society of Photo-Optical 

Instrumentation Engineers (SPIE) Conference Series, Vol. 4836, 
Survey and Other Telescope Technologies and Discoveries, ed. 
J. A. Tyson & S. Wolff, 61-72 
Aldering, C, et al. 2006, ApJ, 650, 510 



20 



Bailey, S., ct al. 2009, A&A, 500, L17 

Baldry, I. K., Glazebrook, K., Brinkmann, J., Ivezic, Z., Lupton, 

R. H., Nichol, R. O, & Szalay, A. S. 2004, ApJ, 600, 681 
Baldwin, J. A., Phillips, M. M., & Terlevich, R. 1981, PASP, 93, 5 
Baltay, C, et al. 2007, PASP, 119, 1278 
Bell, E. F., Mcintosh, D. H., Katz, N., & Weinberg, M. D. 2003, 

ApJ, 585, L117 
Bcrtin, E., & Arnouts, S. 1996, A&AS, 117, 393 
Bloom, J. S., et al. 2012, ApJ, 744, L17 

Bongard, S., Soulez, F., Thiebaut, E., & Pecontal, E. 2011, 

MNRAS, 418, 258 
Borch, A., et al. 2006, A&A, 453, 869 
Brammer, G. B., et al. 2011, ApJ, 739, 24 
Branch, D., & van den Bergh, S. 1993, AJ, 105, 2231 
Bravo, E., & Badencs, C. 2011, MNRAS, 414, 1592 
Bravo, E., Domfngucz, I., Badenes, C, Piersanti, L., & Straniero, 

O. 2010, ApJ, 711, L66 
Brinchmann, J., Chariot, S., White, S. D. M., Tremonti, C, 

Kauffmann, G., Hcckman, T., & Brinkmann, J. 2004, MNRAS, 

351, 1151 
Bruzual, G., & Chariot, S. 2003, MNRAS, 344, 1000 
Bundy, K., et al. 2006, ApJ, 651, 120 
Buton, C, et al. 2012, A&A 
Cardelli, J. A., Clayton, G. C, & Mathis, J. S. 1989, ApJ, 345, 

245 
Chabrier, G. 2003, PASP, 115, 763 
Chotard, N., et al. 2011, A&A, 529, L4+ 
Clcmcns, J. C, Crain, J. A., & Anderson, R. 2004, in Society of 

Photo-Optical Instrumentation Engineers (SPIE) Conference 

Series, Vol. 5492, Ground-based Instrumentation for 

Astronomy, cd. A. F. M. Moorwood & M. lye, 331-340 
Conley, A., et al. 2011, ApJS, 192, 1 

Conroy, C, Gunn, J. E., & White, M. 2009, ApJ, 699, 486 
D'Andrea, C. B., et al. 2011, ApJ, 743, 172 
Davies, R. L., et al. 1997, in Society of Photo-Optical 

Instrumentation Engineers (SPIE) Conference Scries, Vol. 2871, 

Optical Telescopes of Today and Tomorrow, ed. 

A. L. Ardcberg, 1099-1106 
Dilday, B., ct al. 2012, Science, 337, 942 

Domfngucz, I., Hoflich, P., & Straniero, O. 2001, ApJ, 557, 279 
Elbaz, D., et al. 2007, A&A, 468, 33 
Filippcnko, A. V. 1989, PASP, 101, 588 
Filippcnko, A. V., Li, W. D., Treffcrs, R. R., & Modjaz, M. 2001, 

in Astronomical Society of the Pacific Conference Series, Vol. 

246, IAU Colloq. 183: Small Telescope Astronomy on Global 

Scales, ed. B. Paczynski, W.-P. Chen, & C. Lcmme, 121 
Folatelli, G., et al. 2010, AJ, 139, 120 
Foley, R. J. 2012, ApJ, 748, 127 
Foley, R. J., & Kascn, D. 2011, ApJ, 729, 55 
Gallagher, J. S., Garnavich, P. M., Berlind, P., Challis, P., Jha, 

S., & Kirshner, R. P. 2005, ApJ, 634, 210 
Gallagher, J. S., Garnavich, P. M., Caldwell, N., Kirshner, R. P., 

Jha, S. W., Li, W., Ganeshalingam, M., & Filippenko, A. V. 

2008, ApJ, 685, 752 
Gallazzi, A., & Bell, E. F. 2009, ApJS, 185, 253 
Gallazzi, A., Chariot, S., Brinchmann, J., White, S. D. M., & 

Tremonti, C. A. 2005, MNRAS, 362, 41 
Garn, T., & Best, P. N. 2010, MNRAS, 409, 421 
Garnavich, P. M., et al. 1998, ApJ, 509, 74 
Goobar, A. 2008, ApJ, 686, L103 
Gupta, R. R., ct al. 2011, ApJ, 740, 92 
Guy, J., et al. 2007, A&A, 466, 11 
— . 2010, A&A, 523, A7+ 
Hamuy, M., Phillips, M. M., Maza, J., Suntzeff, N. B., Schommcr, 

R. A., & Avilcs, R. 1995, AJ, 109, 1 
Hamuy, M., Phillips, M. M., Suntzeff, N. B., Schommer, R. A., 

Maza, J., & Avilcs, R. 1996, AJ, 112, 2391 
Hamuy, M., et al. 2003, Nature, 424, 651 
Hicken, M., et al. 2009, ApJ, 700, 331 
Hoeflich, P., Wheeler, J. C, & Thielemann, F. K. 1998, ApJ, 495, 

617 
Holtzman, J. A., ct al. 2008, AJ, 136, 2306 
Howell, D. A., Sullivan, M., Conley, A., & Carlbcrg, R. 2007, 

ApJ, 667, L37 
Howell, D. A., et al. 2009, ApJ, 691, 661 



Iben, Jr., I., & Tutukov, A. V. 1984, ApJS, 54, 335 

Jha, S., Riess, A. G., & Kirshner, R. P. 2007, ApJ, 659, 122 

Kasen, D., Ropkc, F. K., & Woosley, S. E. 2009, Nature, 460, 869 

Kauffmann, G., et al. 2003, MNRAS, 341, 33 

Kelly, P. L., Hicken, M., Burke, D. L., Mandcl, K. S., & Kirshner, 

R. P. 2010, ApJ, 715, 743 
Kcsslcr, R., et al. 2009, ApJS, 185, 32 
Kewley, L. J., & Ellison, S. L. 2008, ApJ, 681, 1183 
Kewley, L. J., Groves, B., Kauffmann, G., & Heckman, T. 2006, 

MNRAS, 372, 961 
Knop, R. A., et al. 2003, ApJ, 598, 102 
Kobulnicky, H. A., & Kewley, L. J. 2004, ApJ, 617, 240 

Konishi, K., et al. 2011, ArXiv e-prints 

Kroupa, P. 2001, MNRAS, 322, 231 

Lampeitl, H., et al. 2010, ApJ, 722, 566 

Lantz, B., ct al. 2004, in Society of Photo-Optical 

Instrumentation Engineers (SPIE) Conference Scries, Vol. 5249, 

Optical Design and Engineering, ed. L. Mazuray, P. J. Rogers, 

& R. Wartmann, 146-155 
Lawrence, A., et al. 2007, MNRAS, 379, 1599 
Le Borgne, D., & Rocca-Volmcrangc, B. 2002, A&A, 386, 446 
Lee, J. C, et al. 2009, ApJ, 706, 599 
Lentz, E. J., Baron, E., Branch, D., Hauschildt, P. H., & Nugent, 

P. E. 2000, ApJ, 530, 966 
Li, W., et al. 2011, Nature, 480, 348 
Miller, J. S., & Stone, R. P. S. 1993, Lick Obs. Tech. Rep. 66 

(Santa Cruz: Lick Obs.) 
Ncill, J. D., et al. 2009, ApJ, 707, 1449 
Nugent, P. E., et al. 2011, Nature, 480, 344 
Oke, J. B., et al. 1995, PASP, 107, 375 
Osterbrock, D. E., & Ferland, G. J. 2006, Astrophysics of gaseous 

nebulae and active galactic nuclei, ed. Osterbrock, D. E. & 

Ferland, G. J. 
Pcrlmutter, S., et al. 1997, ApJ, 483, 565 
— . 1999, ApJ, 517, 565 

Pcttini, M., & Pagel, B. E. J. 2004, MNRAS, 348, L59 
Phillips, M. M. 1993, ApJ, 413, L105 
Piro, A. L. 2008, ApJ, 679, 616 
Prieto, J. L., et al. 2007, ArXiv e-prints 
Rana, N. C, & Basu, S. 1992, A&A, 265, 499 
Rau, A., et al. 2009, PASP, 121, 1334 

Riess, A. G., Press, W. H., & Kirshner, R. P. 1995, ApJ, 438, L17 
— . 1996, ApJ, 473, 88 
Riess, A. G., et al. 1998, AJ, 116, 1009 
— . 2004, ApJ, 607, 665 
— . 2007, ApJ, 659, 98 
Salim, S., et al. 2007, ApJS, 173, 267 
Sauer, D. N., ct al. 2008, MNRAS, 391, 1605 
Scalzo, R., et al. 2012, ApJ, 757, 12 
Scalzo, R. A., et al. 2010, ApJ, 713, 1073 
Scannapicco, E., & Bildstcn, L. 2005, ApJ, 629, L85 
Schaefer, B. E., & Pagnotta, A. 2012, Nature, 481, 164 
Schechter, P. 1976, ApJ, 203, 297 
Schlegel, D. J., Finkbcincr, D. P., & Davis, M. 1998, ApJ, 500, 

525 
Skrutskic, M. F., et al. 2006, AJ, 131, 1163 
Smith, M., et al. 2011, ArXiv e-prints 
Spolaor, M., Kobayashi, C, Forbes, D. A., Couch, W. J., & Hau, 

G. K. T. 2010, MNRAS, 408, 272 
Sullivan, M., et al. 2010, MNRAS, 406, 782 
— . 2011, ApJ, 737, 102 
Suzuki, N., et al. 2012, ApJ, 746, 85 
Thomas, R. C, et al. 2011, ApJ, 743, 27 
Timmes, F. X., Brown, E. F., & Truran, J. W. 2003, ApJ, 590, 

L83 
Townslcy, D. M., Jackson, A. P., Calder, A. C, Chamulak, D. A., 

Brown, E. F., & Timmes, F. X. 2009, ApJ, 701, 1582 
Tremonti, C. A., et al. 2004, ApJ, 613, 898 
Tripp, R. 1998, A&A, 331, 815 
Whclan, J., & Iben, Jr., I. 1973, ApJ, 186, 1007 
Wood-Vasey, W. M., Wang, L., & Aldering, G. 2004, ApJ, 616, 

339 
York, D. G., et al. 2000, AJ, 120, 1579 
Zahid, H. J., Kewley, L. J., & Bresolin, F. 2011, ApJ, 730, 137 



